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Heat Storage in Alloy Transformations 


Introduction 

This project has explored the eutectic transformation 
of alloys containing plentiful metals as a means of storing 
heat at moderately high temperatures. In the course of this 
work, the theory of eutectic transformation was examined 
to find guidelines to the best material combinations to 
examine. The heats of transformation were measured 
calorimetrically , and the volume changes of expanding 
solid mixtures and homogeneous liquid solutions, and 
especially during the transformation between the two 
states at fixed temperature, were measured by changes in 
X-ray absorption. Heat flow models appropriate to storage 
in phase change materials have been developed along with 
efficient calculating procedures so that the relative 
importance of the problems associated with energy storage 
density, heat conduction and similar properties can be 
assessed in light of what is possible now. In the plan 
to use alloy eutectics as heat storage media the one step 
that is not completely resolved and demonstrated is 
containment in a manner that will retain the advantages 
of excellent thermal conductivity. It has been demonstrated 
that silicon carbide resists attack by the molten storage 
alloys. However, it has not been demonstrated that carbide 
films can be grown on oxidation resistant alloys and 
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retained through many melting and freezing cycles of 
the storage system. It has been proposed to extend the 
research work in this direction so that working storage 
systems can be designed and tested once containment is 
assured . 

This report is constructed with the cited 
conditions in mind. The first section is a summary 
statement intended to show what has been accomplished in 
preparing heat storage systems utilizing alloy eutectic 
transformations as the working core. The discussion is 
in fairly general terms, because the detailed support 
for the statements is given in a series of extended 
appendices. The appendices are reprints of published 
papers or papers being prepared for publication. In 
several cases more extended discussion of the topics is 
available in Master's or Ph.D. theses accepted by the 
University of Delaware or which soon will be submitted 
to the faculty of the University of Delaware. 

Heat Storage in Eutectic Alloys 

Appendix A, "Heat Storage in Eutectic Alloys," by 
C. Ernest Birchenall and Alan F. Riechman, Met. Trans. A, 
llA , 1415-1420 ( 1980 ) , contains the thermodynamic equations 
governing eutectic transformations in binary systems. The 
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relative costs of elements and their availability indicates 
a current preference for S, Al, Si, Zn, P, Na, Cu, Mg, 

Ca. Fe, although plentiful and cheap, is most likely to 
be needed in the solid containment systems. 

Metal or alloy heat storage by phase change 
appears to be feasible from 1219K to moderately low 
temperatures. Inorganic salts with comparable storage 
capabilities per unit mass or volume also are available, 
as shown in Table II, Appendix A, but the thermal 
conductivities are one to two orders of magnitude higher 
for the metallic storage materials than they are for the 
salts. Simpler and smaller heat exchangers should be 
possible with the metals than with the salts. 

The theory of heats of eutectic transformation shows 
that the important contributions come from heats of fusion 
of the participating elements, possibly modified by their 
combinations in intermetallic compounds, and from the 
differences in the entropies of mixing of the relatively 
well-ordered solid phases and the eutectic liquid. 

Alloys were chosen for study, or sought where they were 
not previously known, in order to optimize the heat storage 
density. For a current tabulation see Appendix B. Although 
it is not possible to state that the best alloys have been 
found, reasonable assurance can be given that the reported 
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alloys are representative of good thermal storage behavior. 
They give the alloys a fair test in comparison with other 
storage materials. 

It is important to consider one other cost factor 
which appears to favor the alloys. Even after 
prolonged use, they remain refined metallic materials 
that have other valuable uses. They retain much of their 
initial value. As a specific example, the Al-Si eutectic 
alloy has been widely used in cast automotive parts, 
including engine blocks. Although many other alloys are 
not in precisely the compositions that would be used in 
other applications, they could be adjusted by simple 
additions. 

Volume Change During Transformation 

To maintain efficient heat transfer from the PCM 
containment vessel to the heat transfer fluid or to a 
heated surface such as the receiver cavity of a solar 
receiver the solid or liquid PCM must always be in direct 
contact. Large volume change during transformation will 
require design complexity to insure that direct conducting 
paths are maintained. The large volume changes reported 
for halide eutectics, on the order of 20 to 25 percent of 
the total volume, were sufficient incentive to measure the 
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changes that occur for the metal eutectics. 


An X-ray absorption method that is readily used 
with ductile samples is described in Appendix C. Its 
application is illustrated by measurements on pure aluminum 
and the Al-Cu and Al-Si eutectics that agree well with 
previously published values. Thermal expansion 
coefficients for the liquid and solid phases or mixtures 
of phases also are in reasonable agreement where published 
values are available for comparison. 

Some refinements have been added to the method in 
Appendix D to simplify the selection of a useful X-ray 
wavelength and improve control of furnace temperature. 

The measurement of fractional volume change during 
eutectic transformation has been extended to two ternary 
systems, with work still progressing. Expansion coefficients 
have been measured again. 

The alloy eutectics show volume changes of about 
5 to 7 pet. for the aluminum alloys and near 3 pet. for 
a Cu-P-Zn eutectic. There is no reason to believe that 
these modest changes are atypical, so the operation of a 
PCM storage system should be simpler with these materials 
than with those that have much larger volume changes 
during transformation. 
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Modeling Heat Plow in Cylindrical Heat Exchangers 
Containing Phase-Change Materials 

Appendix E describes the development of a 
numerical routine for calculating heat flow efficiently 
in a vessel of cylindrical symmetry containing a material 
capable of storing latent heat during a solid to liquid 
transformation, including the special case in which the 
solid is a mixture of phases. This paper was presented 
at the 3rd Miami International Conference on Alternative 
Energy Sources, December 15 to 18, 1980, and will appear 
in the proceedings of that conference to be issued by the 
Hemisphere Publishing Company. The utility of the procedure 
for dealing with metal and salt eutectic storage materials 
and the factors affecting their performances in simple 
heat exchangers is demonstrated in this work. 

An improved procedure is outlined in Appendix P 
that has all of the capabilities of the earlier procedure 
but is capable of more precise calculation. Results have 
been collected with this procedure right up to the 
preparation of the report/ with additional runs to be 
made. The paper is being prepared for the AIAA Thermophysics 
Conference, Palo Alto, Cal., to be held in June, 1981. 

The calculations show the advantages for the highly 
conducting alloys in giving faster storage and retrieval 
of heat for fixed geometry as compared with those materials 
that have lower thermal conductivities. 
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These limited applications have given assurance 
that the calculation techniques can be used in comparing 
alternative designs as well as competitive PCMs in 
establishing optimum system designs. Following demon- 
stration of a proven containment procedure, heat exchanger 
designs consistent with the containment method and storage 
materials will be evaluated, leading to the choice of the 
optimum configurations for the various applications that 
seem reasonable for these combinations. 

Containment 

Recent studies that show the short term stability 
of silicon carbide in contact with molten storage alloys 
are described in Appendix G. Alloys rich in aluminum 
or copper, containing significant concentrations of 
silicon, magnesium and zinc have shown little tendency 
to attack silicon carbide, but an appreciable rate of 
attack on Type 304 stainless steeT (19Cr, 9Ni) . Although 
the exact properties of the containment materials that are 
desirable will be determined by the geometrical con- 
figuration selected, it seems likely that the structure 
will have to be fabricated from an alloy that combines 
good high temperature strength, good thermal conductivity, 
and reasonably good low temperature ductility. In addition, 
resistance to chemical attack by the heat storage materials 
must be provided without impairing these other properties. 
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One reasonable procedure would be to construct the 
system from a stainless steel or other oxidation resistant 
alloy, then coat surfaces that will be exposed to the heat 
storage alloy with silicon carbide. The coating would 
cause little loss in heat transfer and would provide a 
small increase in high temperature strength. If the 
surfaces exposed to air were also coated, somewhat lower 
grade alloys might be used, because silicon carbide should 
provide good oxidation resistance. 

It is important to emphasize that this item is the 
last bit of undemonstrated technology needed to embark on 
a serious design. In particular, the stability of a 
carbide coated steel to many thermal storage and release 
cycles should be verified on a laboratory scale before the 
construction of a working system is undertaken. Several 
substrate alloys, differing in their thermal expansion 
properties, should be coated and tested before a final 
choice of containment material is made. 

A more complete summary of Appendices E and F 
prepared by Michael Labdon follows. 
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Extended Summary of Appendices E and F: Heat 

Transfer with Phase Change Materials in 
a Cylindrical Heat Exchanger 


The thermal storage units depicted in Figures 
1, 2, and 3 represent conceptual designs for use in the 
collection and transmission of energy concentrated from 
parabolic-dish solar receivers. Due to the high energy 
flux imposed on the inner cylinder, the ongoing 
investigation dealt with the contrasting of candidate 
phase-change materials (PCM's) for conducting the heat 
energy, as well as finding an appropriate material for 
PCM containment. 

Due to the array of potential PCM candidates and 
the complicated geometries associated with the storage 
structure, an efficient and versatile numerical routine 
was required for analyzing their thermal behavior. 


Method 

An extensive literature search was initiated to 
determine if a solution method already existed which 
could be modified, or if an original routine needed to 
be developed. The result of the search was to utilize 
a complex method originally posited by Forster (1) , 
Murray (2) and Lazaridis (3) , which was subsequently 
demonstrated as workable by Farouk (4) , and finally 
clarified by Harding (5) and Labdon (6) as to its 
inherent simplicity and efficiency. 
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ih = 1.5 kg/sec 
h = 750 w/m^/°C 



A s 



The versatility of the routine is reflected in the 
fact that time and spatial dependent boundary conditions, 
temperature-dependent properties and initial temperature 
distributions may be included, which represents the more 
realistic condition encountered with solar energy and 
storage systems. Outputs include prediction of melt-line 
locations, temperature histories and charging rates. 

The method models heat conduction within an annulus of 
phase-change material subjected to various boundary 
conditions, by overcoming the inherent difficulty of 
determining the temperature distributions around the 
melt-line (hence, its velocity) which is non-linear in 
nature. Essentially, f inite-dif feirence equations monitor 
the heat transfer in the regions away from the interface, 
while Taylor-series expansions of temperature are made 
near the fusion front so as to describe the temperature 
fields and thus couple the thermal behavior between the 
solid and liquid regions. 

The mathematical derivations, assumptions and 
solution process of the numerical routine appear in two 
publications. Appendix E and Appendix F, where the latter 
models both the radial and axial heat conduction within 
the material. The various thermal properties of the 
candidate PCM's, and the assumptions used in facilitating 
the generation of results are listed on page 15, Appendix E. 
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Analysis of Storage Units 


The solar receiver/storage units dimensioned 
in Figure 1-3 do not represent an attempt in design 
optimization, but merely provide a parity for the contrasting 
of any desirable phase-change materials, which in this 
case are metal alloys and salts. 

Figure 1 was subjected to the simultaneous 
application of heat flux and convection on the inner and 
outer surfaces respectively, with the entire medium 
containing a 10®C initial superheat. In this manner, an 
evaluation of the PCM's can be made while the storage unit 
is in the actual working mode. 

Figure 2 differs from that of Figure 1 by the 
cylindrical array of internal tubes which contain the 
heat-convecting medium. Only the discharging case will 
be modeled by assuming the entire medium is initially ii; 
a constant molten state. Six equally-spaced, 4 cm. diameter 
tubes are situated halfway between the inner and outer 
surfaces. Termination of the discharging case occurs 
when the radially-expanding melt-lines intersect each 
other and the unit's surfaces. The small amount of material 
left over is taken into account by forming an extra concentric 

ring around each tube. 


14 




The third storage unit, depicted in Figure 3, 
represents a more complicated configuration presently 
outside the scope of this work, but it is reasonable to 
expect that the contrasting performance between the salt 
and metal alloys would be similarly evident in this 
geometry as well. 

Results 


Figures 4 through 9 represent output from the 
routine for the first solar receiver/storage unit 
(Figure 1). Figures 4, 5, 6 represent melt-line locations 
within the receiver for the salt, Al-Si, and Mg-Si data 
respectively. The melt-line contours actually represent 
a three-dimensional surface similar to a cone or trumpet- 
end. The more curved lines of the alloys represent the 
increased heat transfer within the upstream nodes due to 
the higher thermal conductivity. The smaller times 
required for melt-line penetration are also indicative of 
the alloys. Interface locations are shown only for the 
outer half of the unit due to the large computer process 
times required for output generation, and due to the fact 
that a steady-state operation is approached where the heat 
flux imposed on the inner cylinder equals the heat 
convected out, thus yielding a melt-line velocity nearing 
zero. 
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Figures 7 and 8 describe the heat flux boundary surface 
temperatures for equal time increments. Both show the 
steady-state temperature distribution being approached, 
as well as the much higher surface temperatures 
associated with the use of the salt PCM; attributable 
again to the low thermal conductivity. Figure 8 indicates 
the effect of increased melt-line penetration on the 
"coQling" or surface temperature reduction of the unit, 
typical of metal alloys. 

A measure of the effectiveness of this type of heat 
exchanger is gleaned from the charging rates shown in 
Figure 9. The area under the curves represents the 
amount of energy transferred within a certain amount of 
time, and as can be seen, the metal alloys are more 
efficient in storing and transferring the collected energy. 

As a further display of the versatility of this 
routine, the discharging rate of Figure 2 was determined 
by the amount of energy transferred from the six tubes, and 
contrasted to that obtained from the discharging case of 
Figure 1 (when the heat flux was terminated) , the results 
of which appear in Figure 10. This indicates the role 
the method plays in design optimization of the receiver/ 
storage units, by showing the increased discharging rates 
associated with installing the convection tubes within 
the PCM, as opposed to a sheet surrounding the entire 
medium.. 
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FIGURE 5 


Melt-line locations at various times (min.) wit 
of the storage unit, for the metal alloy Al-Si. 
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FIGURE 6 


Melt-line locations at eight different times (min.) within the outer 
half of the storage unit, for the metal alloy Mg-Si. 
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FIGURE 7 


Heat flux (inner) surface temperature for the storage unit at 
various times (min.) using the salt phase-change material. 



TIME= 

TIME= 

TIME= 

TIME= 

TIME= 

TIME= 

TIME= 

TIME= 

TIME= 

TIME= 

TIME= 


34.17 

68.17 

102.17 

138.17 

170.17 

204.18 

238.18 

272.18 

306.18 

340.18 
374.1' 


0.60 1.20 1.80 

AXIAL DIS.(CM. ) 



SURFACE TEMP.(C) 

639.00 689,00 739.00 789.00 


FIGURE 8 


Heat flux (inner) surface temperature for the storage unit at 
various times (min.) using the metal alloy, Al-Si. 
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FIGURE 9 


Non-dimensionalized, instantaneous discharging rates of the first 
receiver/storage unit, for the three candidate phase-change materials 
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APPENDIX A 




Heat Storage in Eutectic Alloys 


C. ERNEST BIRCHENALL AND ALAN F. RIECHMAN 

Classical thermodynamic equations based on the regular solution approximation yield 
enthalpic changes for eutectic transformation that agree roughly with values measured for 
several binary and ternary alloy systems by differential scanning calorimetry or differential 
thermal analysis. Restricting measurements to binary and ternary alloys of the relatively 
plentiful elements Al, Cu, Mg, Si and Zn, it has been verified that the best heat storage 
densities on a mass or volume basis are obtained with alloys rich in Si or Al, elements that 
have large heats of fusion. Several of these alloys have the highest heat-storage density 
reported for phase change materials that transform between 780 and 850 K. The MgjSi-Si 
eutectic, which has outstanding storage density at 1219 K, illustrates the utility of ordered 
intermetallic phases with large heat of formation that dissolve in the eutectic liquid to 
contribute to the entropy change. 


Phase change materials (PCMs) are being explored 
as heat storage media to reduce the costs of energy 
generation, conversion and distribution systems. At low 
temperatures, phase changes in Glauber’s salt, Na 2 S 04 
• lOHjO, and other salts have been studied for decades 
as alternatives to heat capacity storage in water or 
rocks.' Rising energy costs and attempts to harvest solar 
radiation have intensified interest in high temperature 
storage. The suitability of PCM systems for this purpose 
has been recognized for years.^ However, systematic 
study of alloys instead of inorganic salts has been 
limited to very rough calculations on pure metals and a 
few binary alloys.^’^ This study explores the character- 
istics that make some alloys better than others. Exper- 
imental heat of transformation measurements on prom- 
ising alloys of the relatively plentiful metals Al, Cu, Mg, 
Si and Zn provide values for the maximum possible 
heat storage densities in these alloys. 

The heat available from a phase transformation 
carried out reversibly at constant temperature and 
pressure is the enthalpy change, which is equal to the 
transformation temperature multiplied by the entropy 
change. A sample calculation applied to the fusion of an 
average close-packed metal illustrates the advantage of 
a PCM over heat capacity storage. 

The entropy of fusion of a normal fee or hep metal is 
shown to be about 10 joules per gram atom per kelvin in 
Fig. !.■* The molar heat capacities are not so uniform, 
but average about 30 joules per gram atom per kelvin. 
Fusion of a gram atom of such a metal at 800 K yields 8 
kJ, which could be stored in heat capacity only by a 
temperature rise of 270 K. In practice, the temperature 
rise is restricted by using much more storage material 
and a much larger container. In many heat storage 
systems, the cost of the containing structure and heat 
exchange surfaces is likely to exceed the cost of storage 
material, so heat storage density is very important. The 
advantages of even larger transformation enthalpies 
that might seem to accrue from the transformation of a 
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condensed phase to a gas are offset by the large system 
volumes that are required to handle the gases. Metals 
and alloys, because of their high thermal conductivities, 
offer an important advantage in that the ratio of heat 
exchanger area to storage volume can be much smaller, 
for a fixed cycling time, than it is for more poorly 
conducting materials.^ 

Any energy storage system having a large central unit 
or many smaller units must use cheap, plentiful mate- 
rials. An analysis^ of a number of elements indicated 
that nine elements, S, Al, Si, Zn, P, Na, Cu, Mg, Ca, 
should be given greatest consideration. Pb addition 
could also be considered in small amounts, while Sb, 

Cd, Sn were ruled out for large scale applications. 

THERMODYNAMIC EQUATIONS 
FOR EUTECTIC MELTING 

Using standard procedures of classical thermody- 
namics and common approximations, equations have 
been written for three cases of eutectic transformation’* 
in binary systems defined by the constructions labeled 
(a), {b), and (c) in Fig. 2. The approximations are: 
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Fig. 1 — Entropies of fusion of some elements. Ga is orthorhombic; In 
is face-centered tetragonal. The figure has been modified from Ref. 4, 
using data from Ref. 8. 
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1) differences between heal capacities of the eutectic 
liquid and solid can be neglected for a moderate 
temperature variation near the eutectic temperature T,, 
or the equivalent assumption, the latent heat L, and the 
entropy change A5 = L/T„ of the eutectic transfor- 
mation are independent of temperature near T/, 2) the 
heat of formation A// from the pure solid elements of 
an ordered intermetallic solid phase t is equal to the 
Gibbs free energy of formation A<7 because the entropy 
change must be very small; 3) the entropy change 
during eutectic melting consists of the entropies of 
fusion of the pure components plus the entropy for 
random mixing of these components in the liquid 
eutectic less the entropies for random mixing in the 
solid solutions. The last assumption is equivalent to 
treating all solutions as regular solutions and all inter- 
mediate phases as completely ordered. All other terms 
in the equation for the Gibbs free energy change can be 
assigned to the enthalpy change, because the Gibbs free 
energy change for the eutectic reaction at 7", must be 
zero. 

Using these approximations, the entropy change for 
case 2 (Fig. 2) is 

AS = — R |(1 — X,) In (1 — X,) -1- X, In x, 

~ In (1 - xj + X, In x„l 

x^^)ln(l - x„) 

+ x, In x,]l -h(l -x„)^+x,^ r,, 

‘A ‘a ^ ' 

Extending this procedure to greater numbers of com- 
ponents adds additional terms of the same kind. In this 
way the entropy change can be increased by mixing 
more and more components if they are present in nearly 
equal proportions. Although this work is concerned 
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Fig. 2 — Eutectic construction for binary systems: (a) Eutectic liquid 
forms from two solid solutions; (6) Eutectic from one solid solution 
and one intermetallic phase; (c) Eutectic from two intermetallic 
phases. 
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with binary systems and a few ternary .systems, it is 
intended to use this base to seek higher order eutectics 
that yield the benefits of increased entropy of mixing. 
The equation also shows that the largest contributions 
are provided by those elements that have the highest 
entropies of fusion. Among the inexpensive elements, 
silicon and aluminum deserve special attention, ac- 
cording to Fig. 1. 

With the entropy contributions identified! the 
enthalpy change at the eutectic temperature can be 
found in two ways: 1) by multiplying A5 by 7',„ ;ind 
2) by evaluating the remaining terms in the free energy 
equation. Comparison of these numbers would be a 
direct check on the validity of the separation of 
contributions if the data were sufficiently precise. 
However, the comparisons also help to identify ques- 
tionable data. The heat effect given by the remaining 
terms in the free energy equation is 

A// = (1 - x,)L, + x,Ls + /?T I (1 - .V,) In r. 

+ >n Yfl - ( VTTt] t 
+ In n] - f -“I-v) t 

+ x^ in [2] 

where the y’s are activity coefficients. 

At the same composition, the activity coefficients 
ordinarily differ more strongly from unity for the solid 
phases than for the liquid phases, a tendency reinforced 
by the increasing deviation from unity with decreasing 
temperature. On the other hand, the most desirable 
liquids are more concentrated solutions than the solid 
phases that react to form them, so the greater deviations 
from unity are likely to be displayed by the activity 
coefficients for the liquids. Unfortunately, the direct 
experimental measurements of thermodynamic quan- 
tities for metallic solutions usually are termin;ited 
considerably above the eutectic temperatures, and often 
are completely unavailable for the solid phases. In those 
cases, only rough extrapolations and judicious guesses 
can be used. Clearly, it is desirable to have large latent 
heats of fusion for all components, combined with 
narrow ranges of solid solubility, in which there are 
strong negative deviations from ideal solution behavior, 
and a eutectic liquid that exhibits much weaker negative 
deviations. 

A second case is found frequently in which a solid 
solution a and an intermediate phase e combine to form 
the eutectic, as in Fig. 2 at point b. Neglecting any 
entropy of mixing in «, the appropriate modifications 
are made readily in Eqs. [I] and [2). 

A5 = - /? I (1 - X,) In (1 - .V,) -f- X, In .v.. 

- [(I - ■'■«) I" (1 - •"■.) + I" '..1 } 

+ (I - + -'"c [-'’1 
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A// = 7? r I (1 - ,v..) In + X, In 

+ (i-,Oi., + (i=-f;)A«, ,4, 

Note that the activity coefficients for the components in 
the intermediate phase are replaced by the equivalent 
heat of formation, which is more likely to be tabulated. 

The final case to be considered is that in which two 
ordered intermediate phases e and to combine to form a 
eutectic melt as illustrated in Fig. 2 at point c. The 
analogs of Eqs. [1] and [2] now become: 

AS = -/? [(1 - x,)ln(l - X,) 

+ X, In x,l + (1 - X,) -^ + X, ^ 

IH = RT [(I - X,) In Y' + x. In yi] + (1 - x,)L^ 

[ 6 ] 


The last case makes it clear by comparing Eqs. [5] 
and [6] that the separation of entropy and enthalpy of 
transformation cannot be exact, because Eq. [5] does 
not contain any explicit properties of t and to. They 
only need to form with sufficient strength of binding to 
suppress substitutional disorder. Thus, Eq. [6} seems 
likely to give a better estimate of the heat storage 
capabilities of such systems. In the calculations to 
follow, the values for estimated heat storage based on 
TAS, that is on Eqs. [1], [3] and [5], are designated as 
the first approximation, while the values based on A//, 
as given in Eqs. [2], [4] and [6] are designated as second 
approximations and are given in parentheses in Table I 
The latter, however, are more sensitive to experimental 
errors in the tabulated data, so they are not necessarily 
more reliable than the first approximations. Our direct 
measurements can be compared with these calculated 
values. 

Experimental Procedure 

The enthalpy changes for several binary and ternary 
eutectic alloys were measured by differential scanning 
calorimetry for temperatures below 1020 K (Perkin- 
Elmcr DSC-2) and by differential thermal analysis for 
temperatures above 1020 K (Du Pont 900 DTA with 
high temperature cell). The compositions and temper- 
atures of transformation were obtained from published 
phase diagrams.’"'® Alloys were melted from commer- 


Tabla I. Heats at Eutectic Transformation (Maximum Heat Storage Capacity) and Undercooling ol Selected Binary and Ternary Metal 

Eutectics 


Eutectic Alloy, 
•Mol Fractions 

Eutectic 

Undercooling, 

K 

Compositions’ 

Heat of Formation,” 
kJ/g • at. 

Maximum Heat Storage, ki/kg 

Temperature, K 


Calculatedt 

Measured 

Binaries 
Mg-0.29 Zn 

616 

2-4 

0.025 0.30* 

12.6 

247 

138 






(464) 


AI-0.375 Mg 

724 

5-9 

0.19 0.38J 

3.05 

458 

310 






(477) 


AI-0.175CU 

821 

1-3 

0.025 0.333 

13.3 

359 

351 






(380) 


Al-0.13Si 

852 

2-3 



571 

515 

Mg-0.529 Si 

1219 

1-2 

1.0 0.333 

26.4 

1212 

774 




Proportion"* 
of Phases 




Ternaries 

A1-0.17CU 

779 

4-10 

/.= 0.354 

r=;L5 

406 

360 

-0.162 Mg 



( 8.0 
V 13.25 

(400) 





/oiMraii “ 0.620 






/coAii 0.026 

=j 3.6 (eutectic) 



AI-0.126 Cu 

833 

5-6 

/mr,s. = 0.067 

26.4 

549 

545 

-0.051 Mg 




(449) 





fs, = 0.094 





f 0.981 Al 

/. = 0.838 <0.11 Si 

1,0.0076 Mg 


• Values not in parentheses are calculated from Eqs. [1], [3], or [5). Values in 
parentheses are calculated from Eqs. [2], [4] or (6j. 

■f Structure of intermediate phase is unknown. Heat of formation estimated to be less 
than that of MgZn^flV.b kJ/g • at.). An even smaller A/f, or a disordered solid phase 
would lower the calculated heat storage in parentheses. 

I Melts congruently at AljMg; and 724 K. 
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cially pure elements in graphite crucibles. The eutectic 
structure, in each case, was verified by optical metallo- 
graphy to be fine enough that a small sample would 
represent the bulk alloy. Proper behavior of the DSC 
and DTA curves also verified the absence of any 
significant excess of either phase. 

The measurements were calibrated by runs on six 
carefully purified elements. Small alloy samples were 
melted and refrozen in graphite liners to insure good 
thermal contact prior to the first measurement. The 
samples were then measured repeatedly to ascertain 
that the transformation temperatures and heats of 
transformation were reproducible. For the DSC meas- 
urements, good determinations of heat capacity also 
were made. Only the MgjSi-Si alloys at very high 
temperature in the DTA apparatus showed a 'drift in 
behavior, which was attributed to small losses of the 
volatile Mg metal. Although this drift was minimized by 
using crucibles with good closure in inert gas, the values 
reported are those obtained in the first cycle. Even so, 
the enthalpy change listed for that alloy may be low. 

The alloys selected for measurement are those con- 
taining the plentiful elements for which the theory and 
preliminary calculations based on published data pre- 
dicted the highest heats of transformation. Other eutec- 
tics from these same systems and eutectics from other 
binary systems did not seem likely to yield maximum 
storage densities as high as these. Many ternary systems 
are incompletely known, so other ternary alloys may be 
measured later from ternary systems that have been 
examined in part. 

Calibration of the DTA apparatus indicated that a 
precision of ±4 pet and ± 3 K could be expected for 
the enthalpy of fusion and temperature, respectively. 
Heat capacity could not be obtained with useful pre- 
cision. For the DSC apparatus, precision of ±2 pet and 
± 1 K could be expected for the enthalpy change and 
temperature, respectively. Furthermore, heat capacities 
could be measured to about ± 10 pet or better. Addi- 
tional uncertainties can arise if the composition is not 
precisely at the eutectic value, but a small amount of an 
excess phase makes a comparable heat contribution of 
its own, so that the effect should not be a .sensitive one. 

Summary of Calculations and Measurements 
on Binary and Ternary Eutectic Alloys 

Table I summarizes measurements and calculations 
for the chosen systems. For Mg-Zn and Mg-Si, the 
measured values lie far below the calculated values. The 
Al-Mg result also is seriously below predictions. Al-Cu 
and Al-Si-Mg give very close agreement. Al-Si and 
Al-Cu-Mg are satisfactorily close to the predicted 
values. Calculations based on the available, tabulated 
data are not a satisfactory substitute for measurements 
on the eutectic alloys. 

The undercoolings recorded for all of these alloys are 
based on the very small samples used in the thermal 
analysis measurements. The starting materials were 
commercially pure metals, not research grade materials. 
However, commercially pure metals usually contain 
very little impurity to act as nucleation catalysts. The 
small samples tend to undercool more than larger 
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amounts of the same alloy. If there is a need to restrict 
undercooling to less than the 5 to 10 K range observed 
for several systems, intentionally added dopants and 
larger samples should eliminate the problem. In any 
case, a real system is not likely to undercool to the 
extent reported here. 

The heat capacities measured on the DSC for all 
alloys listed in Table I are given in Table II. Those heat 
capacities apply to the solid mixture at the transfor- 
mation temperature. Some of the other heat capacities 
in Table II taken from previous tabulations may apply 
at room temperature. 


DISCUSSION 

Definitive comparisons among the materials that 
have been proposed for high temperature storage are 
premature because the tabulated data may be uncertain, 
and some critical data are completely missing. Table II 
contains a preliminary comparison with salts identified 
by Tye et a!,'' and Kauffman and Lorsch^ to have merit 
because of unusually high storage density or a good 
temperature match to a particular heat generating 
device. Lithium-containing compounds and fluorides 
stand out because of their low atomic numbers. The 
studies of Schroder'^ have emphasized fluoride eutec- 
tics, and the Dynatech studies" have emphasized 
lithium compounds and eutectics. 

The relatively high cost of these materials, especially 
of lithium and lithium hydride, must be borne in mind. 
Also, the salts must be dried very carefully and kept drv 
to limit their corrosive attack on containers. Volume 
change measurements on transforming salts and salt 
eutectics are sparse. However, their volume changes for 
solid-liquid transformations sometimes are large;^ for 
example, values in excess of 20 pet are reported for 
melting LiF, LiCl, NaF, NaCl, LiNOj, and the 
63LiOH/37LiCl eutectic." In contrast, the alloy trans- 
formations appear to involve changes appreciably less 
than 10 pet. For the AI-AljCu eutectic transformation 
the volume was found to increase by 5.1 pet during 
fusion by Birchenall, Harrison and Balart.'^ 

Where thermal conductivities of salts are given, they 
lie in the range 1 to 5 W/m • K, except that LiH has a 
value of 7. The solid metals have thermal conductivities 
of 40 to 400 W/m • K, and the alloys and liquid metals 
may drop to about half of these values, retaining an 
advantage of one to two orders of magnitude. 

Tyec/a/" specified important temperature ranges to 
match the principal heat-generating sources: water- 
cooled nuclear reactors (PWR and BWR), 505 to 
544 K; fossil-fueled supercritical-steam reactor (FFR), 
783 to 853 K: high-temperature, gas-cooled, graphite- 
moderated nuclear reactor (HTGR), 977 to 1033 K. The 
storage materials studied for these cases were: LiNOj 
(527 K for PWR), 63LiOH/37LiCl eutectic (533 K for 
BWR), LiOH (743 K for FFR), and NajB,0, (1015 K 
for HTGR). However, Na 2 B 40 ;, usually did not crys- 
tallize. Glass-forming borates, silicates and phosphates 
are unsuitable PCM’s for heat storage applications. 

Thus, lithium salts are the basis of all of the workable 
materials in that study. 
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Tabit II. Comparison ol l>hase-Changa Materials as Heat Storage Media 


Material 
(wi pci) 

Transformation 

Temperature, 

K 

Heat of 

Transformation, 

kJ/kg 

Solidt 

Density, 

kg/m’ 

Volumetric^ 
Heat Storage, 
kJ/m’ 

Solid Heat 
Capacity 
kJ/kgK 

Coefficient of* 
Thermal 
Conduction 
W/mK 

Mg.Si/Si 

1219* 

774* 

(2000)t 

1.55 X 10“ 



26/20 

LiF 

1121 

1044 

2640 

2.76 

1.64 

1.92 

75NaF/25MgF 

1105 

650 

2680 

1.74 

1.42 

4.66 

67LiF/33MgFj 

1019 

947 

2630 

2.49 

1.42 

. — 

65NaF/23CaFj/12MgFj 

1018 

574 

2760 

1.58 

1.17 

— 

Li,CO, 

998 

605 

2200 

1.33 

2.64 

— 

CaMg, 

990 

554 

(1550)t 

0.86 

— 

4/20 

MgClj 

988 

454 

2240 

1.02 

0.75 

— 

LiH 

956 

2582 

790 

2.04 

8.04 

7.0 

A1 

934 

400 

2370J 

0.95 

1.29* 

204.2 

33.4LiF/49.5NaF/ 17.1 MgFj 

923 

860 

2810 

2.42 

1.42 

1.15 

Mg 

922 

368 

1590t 

0.59 

1.34 

131 

46LiF/44NaF/10MgFj 

905 

858 

2610 

2.24 

1.40 

1.20 

MgZiij 

861 

275 

5200 

1.43 

— 

— 

Al/Si 

852‘ 

515* 

(2250)* 

1.16 

1.49* 

180/70 

Ca(NO ,)2 

834 

130 

2500 

0.33 

0.88 



Al/Si/Mg 

833* 

545* 

(2300)* 

1.25 

1.39* 

200/70 

AI/A1;CU 

821* 

351* 

3424t 

1.20 

1.11* 

130/80 

Mg/Mg,Ca 

790 

(353) 

(1570)t 

0.55 

— 

_ 

Al/AljCu/AljCiiMg 

779* 

360* 

(3050)t 

1.10 

1.09* 

115/75 

5f)Na.CO/44Li;CO) 

769 

368 

2330 

0.86 

1.85 

2.1! 

LiOH 

744 

1100 

1340t 

1.47 

4.5 

1.3 

Al/MgjAl, 

724* 

310* 

(2300)t 

0.71 

1.73* 

80/50 

50NaCl/50MgClj 

723 

429 

2240 

0.96 

0.93 

0.96 

LiOH/LiF 

700 

1163 

1150 

1.34 

2.14 

1.2 

Zn 

693 

112 

7140 

0.80 

0.45 

— 

3 1 LijC03/35KjC0/33N,ijC0j 

670 

275 

2310 

0.64 

1.69 

2.04 

63MgCl,/22.3NaCl/ 14KC1 

658 

461 

2250 

1.04 

0.96 

0.95 

.Mg/Mg,Zn 

616* 

138* 

(4900) 

0.68 

1.04* 

80/50 

NaOH 

593 

160 

2070 

0.33 

1.47 

1.54 

95.4NaNOj/4.6NaCl 

570 

191 

2260 

0.43 

1.85 

0.61 

7.*NaCI/6.4NajCO,,''85.8NaOH 

555 

316 

2130 

0.67 

2.51 

— 

37LiCl/63 LiOH 

535 

485 

1550 

0.75 

2.4 

1.10 

LiNOj 

527 

530 

2120 

1.12 

2.05 

1.37 

Li 

452 

663 

534 

0.35 

3.98 

— 


• Measured or confirmed in this work. 

t Density values in parentheses are estimated from the behavior of the elements. Other values are from tables cited, 
t Density at transformation temperature; others appear to be at room temperature, which inflates volumetric heat storage density. 
■‘These values depend on the densities. Nearly all will be lowered further when based on liquid volume. 

t Where two numbers are given, Che first is an estimate for the solid, the second for the liquid, both at transformation temperature. 


Other large-scale applications may be possible in 
these and other temperature ranges. For example, 
central station solar thermal power generation might be 
compatible with storage between 600 and 930 K, 
depending on the type of conversion cycle with which it 
is coupled. Intermediate scale industrial and transpor- 
tation applications also may spread over a wide tem- 
perature range, including temperatures above 1600 K 
for steel soaking pits. Small-scale storage applications, 
at widely scattered temperatures, also take place down- 
stream from the power generator. If those applications 
are built around a heat-storage system, they can level 
the load on the generation and distribution systems. For 
that reason, Table 11 is not limited to the three ranges 
emphasized by Tye et «/." Entries are listed according 
to transformation temperature. Values enclosed by 
parentheses are our calculated or estimated values. 
Those not so enclosed are either measured values, or are 
taken directly from tables.^'" The density values with 
double daggers are for the solid at the transformation 
temperature. All other densities appear to be room tein- 
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perature values, which gives corresponding volumetric 
heat storage values that are slightly too high. The solid 
heat capacities, as well as heats of transformation, and 
transformation temperatures measured by DSC in this 
investigation are distinguished by asterisks. 

CONCLUSIONS 

From Table II, it appears that the Al-Cu, Al-Si, 
Al-Cu-Mg and Al-Si-Mg eutectics have the best heat 
storage characteristics in the range to operate with 
fossil-fueled combustors. The MgjSi-Si eutectic has very 
good heat storage capacity at 1219 K, where it might 
find application in a solar power tower. Heat conduc- 
tion should offer added advantages over other types of 
materials, especially for short charging or discharging 
times.^ In the HTGR and central-station, solar-thermal 
range, the metallic materials must compete with fluoride 
mixtures unless the operating temperature is high 
enough to be compatible with the magnesium-silicon 
eutectic, or unless that transformation temperature can 
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be lowered by further alloying. In the intermediate 
temperature range LiH has outstanding storage capa- 
bility on a mass basis, but it is very reactive and 
expensive. Ternary fluoride mixtures appear to be 
preferable on a volume basis. 

A similar sort of comparison might be made between 
LiOH, one of the least expensive lithium compounds, 
NaCl-MgClj, and the AI-Mg eutectic, where high ther- 
mal conductivity may favor the alloy for some purposes. 
LiOH/LiF, LijCOj/KjCOj/Na^COj, and Zn may offer 
a similar choice, to be determined by conduction and 
containment problems that are not definable without 
specifying an application. If metals are to play any role 
below the Mg/MgjZn eutectic temperature, it must be 
for small-scale special applications that justify the use 
of relatively expensive metals, or they must wait for 
ternary or more complex systems that remain to be 
identified. 
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ABSTRACT 


Eutectic alloys and interraetallic compounds of the 
common elements, Al, Ca, Cu, Mg, P, Si, Zn, may be useful 
for high temperature heat storage. A number of binary and 
ternary eutectic alloys have already been studied and their 
heats of fusion measured. In this work heats of fusion, and 
in come cases heat capacities, of new multicomponent 
eutectics and intermetallic phases are reported. Some of 
the alloys investigated are competitive with molten salts 
for heat storage density at high temperatures. 

A method to determine unknown eutectic compositions 
was developed, based on combined results of differential 
thermal analysis with metallography and microprobe analysis. 
A new differential thermal analyzer, which gave satisfactory 
accuracy for qualitative measurements, was built to handle 
the relatively large samples that are required for this 
procedure. The method allows determination of eutectic 
composition in no more than three steps. 

The heats of fusion of the resulting alloys were 
measured using commercial calorimeters, a differential 
thermal analyzer and a differential scanning calorimeter. 

The possibility of measuring heat capacities by 
differential thermal analysis was also examined. 



INTRODUCTION 


There are various ways in which heat storage can 
reduce costs and make more efficient use of power 
generating equipment. A thermal storage unit would 
eliminate the difference between peak and off-peak hours, 
storing heat when demand is less than the generating 
capacity and returning it as needed. The generator could 
be designed for the average, instead of the maximum, power 
requirement. A central storage system could thus reduce 
the cost of building and operating a power generating plant. 

An efficient way to store heat is to use phase 
transformations that take place at a definite temperature. 
The storage density depends upon the latent heat of 
transformation, which is the entropy change times the 
transformation temperature. This entropy change is small 
for solid-solid transformations but much larger in melting. 
For solid-gas transformations, advantages of the still 
larger entropy change are offset by the practical problems 
in handling the large volume of gases. 

The first materials considered for high temperature 
heat storage were inorganic salts. Fluoride eutectics 
were studied by Schroeder (1) . Lithium compounds and 
eutectics were characterized by Tye et al. (2). Further 
study in molten salts was carried out by Ferrara et al. (3), 
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by Haru et al . (4) , and by Clayton (5) . Some of the 

problems found in molten salts were: 

i) They have low thermal conductivity, leading 
to large irreversible losses and requiring 
a very large surface to volume ratio. 

ii) They are highly corrosive and must be dried 

very carefully, even if expensive, corrosion- 
resistant containers are used. 

iii) They have large volume changes during melting, 
around 20%. 

iv) Lithium and fluoride compounds, those with the 
best heat storage densities, are relatively 
expensive materials. 

Metals diminish these problems, and therefore might 
be strongly competitive even if they have somewhat lower 
storage densities. Materials were selected in the range 
of 300 to 1000°C, on the basis of cost and availability, 
as well as storage density (6). The heats of fusion of 
many binary and ternary eutectics were calculated 
approximately (7) , and the most promising of these 
were measured (8) . 

In choosing particular eutectics, calculations were 
done for the heats of fusion using two different ideal 
solution models as well as a regular solution model. 
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These calculations are discussed below. The results 
are only approximate, but it was possible to identify 
qualitatively the contributions to the heat of fusion 
in order of importance. 

i) Heats of fusion of the constituent elements, 

ii) Heat of formation of intermetallic compounds 
present. 

iii) Difference in entropies of mixing between 
solid and liquid phases. 


Analysis of the results obtained by Riechman shows that 
alloys have the best observed heat storage densities in 
the range from 700°K to 850°K. Better alloys are sought 
in the region 900°K to 1000°K and below 700°K by further 
alloying. Therefore, it is necessary to determine the 
congruently melting compositions in multicomponent systems 
and evaluate the heats of fusion. 


Ternary phase diagrams are not generally known, 
and few thermal data are available. In some cases 
the formation of ternary intermetallic compounds is reported 
With this lack of information any theoretical approach to 
the problem is a crude approximation. Nevertheless, such 
an approach can be the useful starting point of a trial- 
and-error method for experimental determination. 



Although it is clear that alloys will be used primarily 
as latent heat storage systems and only incidentally for 
sensible storage, it is important to know their heat 
capacities as well. Complete information is needed to 
analyze the possible heat storage systems and develop 
practical design. 
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CALCULATION OF EUTECTIC COMPOSITIONS AND TEMPERATURES 

For the ideal solution model all calculations are 
simple, and no thermal data are needed except heats of 
fusion of the Components. For the regular solution case, 
the binary equations have been extended to ternary systems, 
an approximation that allows calculation of ternary free 
energies from binary data. This was first done by Darken 
(9) in 1950 and developed further by Toop (10) . The 
derivation is very well presented by Richardson (11) . 

An alternative but equivalent treatment is given by 
Olson and Toop (12) . This approximation has been used 
successfully in the calculation of phase diagrams with no 
ternary intermetallic compounds, for example, the Ni-Cu-Cr 
system by Meijering (13) . The limitations of the method 
were analyzed by Olson and Toop (14) . In cases where 
compound formation is suspected, they concluded that the 
regular solution model does not yield useful results. 

In addition to these two equations, a polynomial 
approach has been used. This Guggenheim-type equation has 
been generalized by Hillert (15) . 

The three approaches were compared by Ansara et al. 

(16) for the determination of phase boundaries. In the 
systems Fe~Cr-V and Cd-Bi-Sb it was found that the results 
derived from these three equations did not differ much. In 
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cases where compounds are formed or where there is a strong 
ordering in the liquid these approaches are not expected 
to be useful, especially the first two models that assume 
regular solution throughout the system. The advantage of 
these models is that they are based on binary data and no 
ternary measurements are needed. 

In this work both ideal and regular solution models 
were applied and the results compared. Clausius-Clapeyron 
equations, one for each component phase, and the condition 
that all the fractional compositions sum to unity yields 
n + 1 equations in n unknown compositions and the 
temperature. This set of non-linear equations was solved 
by computer fitting. The predicted eutectic temperatures 
for some alloy systems are shown in Table 1. 

A regular solution model applied on the basis of 
pure components will not give adequate results if 
intermetallic compounds are present in the system. For 
this reason the same approximation was used taking the 
constituent solid phases as the components of the system. 
Thermal data for the calculations were taken from 
Hultgren et al. (17), except in the Mg-Si case (18), 
where activities are always given with the pure elements 
as standard states. These reference states were changed 
to be the pure intermediate phases. The eutectic points 
predicted by this model are also shown in Table 1. 
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Table 1, Calculation of Eutectic Compositions 
and Temperatures 


System 


Ideal 

Model 


Regular 

Model 


Experimental 

Al-Cu 

15% 

Cu 

654 

°K 

17% 

Cu 

790 

°K 

17% 

Cu 

821°K 

Al-Mg 

36% 

Mg 

555 

°K 

40% 

Mg 

676 

°K 

37% 

Mg 

724“K 

Al-Si 

1 

y% 

Si 

— 


7% 

Si 

834 

°K 

12% 

Si 

850°K 

1 

Mg-Si 

52% 

Mg 

1183 

°K 

53% 

Mg 

1289 

°K 

48% 

Mg 

1219°K 

Mg-Zn 

32% 

Zn 

592 

°K 

33% 

Zn 

655 

°K 

29% 

Zn 

615°K 

Mg-Al-Cu 

10% 

Mg 18% 

Cu 546 

°K 

11% 

Mg 17% 

Cu 541 

°K 

8% 

Mg 

17%Mg 797°K 





CALCULATION OF THE HEATS OF FUSION 


For eutectic melting the ideal solution model can 
be applied to the prediction of heats of fusion using either 
the participating elements or the constituent phases as 
components. The two approaches yield different results. 

Both assume complete ordering in the solid with different 
ordering in the liquid. Random mixing of elements is 
likely to be better for disordered liquids and random 
mixing of phases for partially ordered liquids. 

Regular solution equations for eutectics (19) 
permit calculation of the latent heat either from the 
entropy change or from the enthalpy change, assuming 
random mixing in the liquid. If the solutions are truly 
regular and the data are correct, both equations should 
be equivalent. 

Table 2 is a comparison between the described 
theoretical methods and experimental results. It 
appears that the ideal mixing calculation from the parti- 
cipating phases often is the best choice for accuracy and 
simplicity. The calculations by all methods tend to 
overestimate the heat of fusion, and they should be taken 
as an upper bound to the true value. 
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Table Z. Calculation of Heats of Fusion 




AH^ (cal/gr) 

Alloy 

T(°K) 

Ideal 
1 * 

Ideal 
2 * 

Regular 
1 * 

Regular 
2 * 

Measured 

Al-Cu 

821 

106.3 


86 

91 

83.8 

Al-Si 

852 

136.5 

136.5 



123.1 

Al-Cu-Si 

797 

134.0 

97.8 



101 

Al-Mg-Si 

833 


137.1 

131 

107 

130.3 

Al-Mg-Zn 

720 


95.9 



74 

Al-Mg 

724 


89.6 

109 

114 

74.1 

Mg-Cu-Zn 

725 

97.5 




60.5 

Mg-Zn 

615 


55 

59 

111 

40 

Al-Mg-Cu 

779 


89.7 

97 

96 

86.1 

Mg-Ca 

790 

102,8 




63.1 

Mq-Si 

1253 

300 

256 



185 


* In column 1 the pure elements are taken as components; in column 2 the 
participating phases are taken as components. 



EXPERIMENTAL DETERMINATION OF EUTECTIC COMPOSITIONS 

An experimental method for finding and characterizing 
eutectics was needed that would work with no information at 
all on the system but would be able to use known information 
to converge faster. The combined results of DTA, electron 
microprobe and metallography provide such a trial-and-error 
method. Any theoretical approximation or empirical data 
can be used to improve the initial guess of eutectic compo- 
sition. 

A sample of an arbitrary or guessed starting composition 
is prepared, and the melting behavior is determined by DTA. 
Metallographic inspection of the alloy is done optically 
and in a scanning electron microscope (SEM) . The observed 
phases are correlated with the DTA results. The sample is 
searched for a small eutectic region, and the composition of 
any small zone suspected of being a eutectic mixture is 
analyzed quantitatively by SEM with an energy dispersive 
analyzer for the emitted x-rays (EDAX) . If some doubt 
remains about the number or identity of the phases present, 
x-ray diffraction measurements are made. 

If there is a eutectic region with a melting tempera- 
ture in a desirable range, a new sample is prepared with 
the composition obtained by EDAX on the previous sample. 

The procedure is repeated until an alloy that is completely 
eutectic is obtained. When an intermetallic compound is 
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found in the process, it is also studied, especially in 
higher order systems where the probability of finding a 
desirable eutectic is lower. 

In order to be sure that the observed DTA peaks corres- 
pond to the observed phases and the compositions measured, 
it is necessary to work with the three techniques on the 
same seimple. In this case inhomogeneities in the alloy are 
not a serious problem. 

Commercial calorimeters v/ork with very small seimples in 
order to minimize the errors in quantitative measurements. 

A DTA apparatus that will give good qualitative measure- 
ments with a large sample is diagrammed in Fig. 1. Temper- 
ature was varied under control of a programmer. The 
differential temperature was amplified and plotted directly 
against sample temperature on an XY recorder. 
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PREPARATION OF ALLOYS 


The first alloys were made in a resistance furnace 
following the procedures of Riechman (8,19) . However, 
when the starting materials differed greatly in melting 
temperatures, vapor losses could not be avoided, even with 
a collar on the stirrer shaft closely fitted to the crucible. 

Owing to the repeated failure to prepare the 
alloys in which the starting materials have large 
differences in melting temperatures, an induction furnace 
was assembled. The relatively low frequency (50 kHz) of 
the mercury spark gap resonator provided vigorous internal 
stirring to the melt. 

The starting materials were introduced into alumina 
crucibles that were then sealed by an alumina disc of 
the same diameter cemented with Sauereisen cement No. 78. 

In the cases where most of the starting materials were semi- 
conductors, graphite crucibles were used, sealing them in 
the same way. The crucible was placed in a mullite tube, 
with a thermocouple touching it from the bottom and argon 
flowing in from the top. This setup was entirely successful 
in preventing vapor losses and providing good stirring to 
the melt. 
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MEASUREMENT OF HEATS OF FUSION AND 
HEAT CAPACITIES 

Two commercial instruments were used for the 
measurement of heats of fusion, the Perkin Elmer DSC-II 
and the Du Pont 990 DTA. The former gave better precision, 
about 3%, and was used for alloys with melting points 
below 700°C. The Du Pont 990 DTA with its high temperature 
cell gave measurements with 5% accuracy for temperatures 
above 700°C . 

For alloys with melting temperatures below 700°C 
the measurement of heat capacities can be carried out by 
DSC with a 10% accuracy. Manufacturers of commercial 
DTA instruments for measurements at higher temperatures 
do not claim sufficient baseline control for heat 
capacitv determination. However, the Du Pont 990 DTA was 
tested for measuring heat capacities because of its 
improved electronic control of the baseline. Twelve 
baselines were analyzed, and it was found that they were 
straight within 10%. To check reproducibility several 
runs were made with the same sample pan that was opened 
and repositioned in the furnace between runs. 

After two successful calibrations with sapphire, 
measurements on Al and Ag gave results shown in Table 3. 
For Al the agreement with tabulated data is very good. 

For Ag agreement is not as good, because its heat capacity 
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Table 3 


Heat Capacity of Al and Ag 
Measured with DTA 1200°C Cell 


T(C°) 

Cp (meas) 

Cp (tabulated) 


(cal/g) 

(cal/g) 

Aluminum 

550 

.26 

.28 

580 

.26 

.28 

600 

.26 

.28 

630 

.27 

.29 

650 

.29 

.29 

(liq) 710 

.36 


Silver 

800 

.12 

.07 

850 

.10 

.07 

900 

.15 

.07 

930 

.25 

.07 

(liq) 1000 

.5 

- - 


B16 














taken at higher temperatures, is much lower. At lower 
temperatures the agreement improves. 


Measurements were also attempted for the Mg-Si, 
Cu-Si and Mg-Cu-Si eutectics. At high temperatures, the 
values obtained are too low. In some cases the solid and 
liquid baselines crossed, owing to the evaporation of 
the volatile elements. 

It appears that the specific heat values from DTA 
for non-volatile alloys are reliable within about 20%, 
but the values for volatile alloys cannot be trusted. 
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RESULTS 


New Eutectics and Intermetallic Compounds 

A search for new eutectics and intermetallic 
compounds was conducted in several ternary systems and 
two quaternaries. The search method converges in three 
steps, in two cases to negative results. The diagram 
in Table 4 summarizes the results step by step. 

X-ray data on the new compounds found are reported 
separately (20) , together with several interesting 
microstructures . 

Measurement of Heats of Transformation 

Heats of eutectic transformation were measured in 
the Perkin Elmer DCS II when the melting temperature was 
lower than 700 °C. For higher temperatures the Du Pont 
DTA model 990 was used with its high temperature cell. 

Several known eutectics were characterized. 

Two of the previously measured eutectics (19) were 
repeated, Mg-Zn because the previous result was believed 
to be too low, and Mg-Si for further check. For Mg-Zn a 
higher value was obtained, and for Mg-Si excellent 
agreement was found with the previous measurement. The 
results for the known alloys are listed in Table 5, 
together with those newly determined. Table 5 also 
includes the measurements of heat capacities for the 
alloys that are interesting as heat storage materials. 
The room temperature densities used in calculating 
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Table 4 . 


StoiKi in tbo Search for N'ew Eutectics and Intermetallic Compounds 


£ys ter. 

' 

Alloy No. 1 

M 1 c r o 3 1 r u c t. u r t': Ti . e r mog r am 

Alloy Mo. 2 

Mi orostructure Thermogram 

Alloy No. 3 

Microstructure Thermogram 

Mg-3i-Cu 

.’■••ggSigCu^g^ 

segregated 

2 peaks 

eutectic 

2 peaks 

ternary 

eutectic: 

«ggSi7CUig, 

Mg2Si,Si 

1 peak 
770*C 

Hg-Cu-Ca 

ternary 
eutectic + 
segregated 
Mg^Cu 

1 peak 
453*, bump 
after the 

ternary 
eutectic: 
Mg^Ca, Mg, 

ternary comp. 

1 peak 

445*C 



Mg-Si-Zn 

MgSi eut.+ 
an 31 Mg 
eutectic 

1 broad 
peak 
800 *C 

Mg-Si eut. 
+Zn 3% eut. 

1 broad 
peak 
850*C 

Mg-Si eut, 
+ Zn 

1 broad 
peak 
900*C 

P-Cu“Zn 

eutectic : 

6 brass, 

CUjP, 

ternary comp. 

1 peak 
720*C 





P-Cu-Si 

1 hcmogeneouD 
phase 

1 peak 
840*C 





Cu-Si-Zn 

ternary 
compound + 
ternary 
eutectic 

2 peaks 

ternary 

eutectic: 

Cu^Si, 

ternary comp, 
otlier phase? 

1 peak 
765*C 




CaSi^ + Si 
+ Ca2HgSij 

2 peaks 

CaSi^ + Si 
+ Ca2MgSi^ 

2 peaks 

Mg2Si + Si 
+ ca2,MgSi2 

1 peak 
• 65*C 

Al-C u-Mg 
-Si 

Al^Cu t Si 
+ Mg., Si + 
quat. comp. 

4 peaks 

AI2CU + Si 
+ Mg, Si t 
quat comp. 

'4 peaks 



Al-Cu-Mg 

2 quat . comp. 
+ quat. 
eutectic 

2 peaks 

eutectic; 

2 quat .comp., 
Al, other 
phase 

1 peak 
520*C 



Al-Mcj-Si 

-2n 

Al-Mg-Zn eut. 
+ Mg^Si + 

Si 

2 peaks 

Al-Mg-Zn eut. 
+Mg2Si+Si 

2 peaks 
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Table 5. Summary of Heat of Fusion Measurements 


System 


Mg-Zn 

MgZn2 

Mg-Si* 

Cu-Si* 

Mg-Ca 

Mg^Cu 

Cu-P 

Al-Cu-Si* 

Mg-Cu-Zn LT 

Mg-Cu-Zn HT* 

Al-Mg-Zn* 

Ca-Mg-Zn 

Cu-Mg-Si* 

Mg-Cu-Ca 

Mg-Si-Zn HT 

P-Cu-Zn* 

P-Cu-Si 

Cu-Si-Zn 

Mg-Si-Ca 

Al-Cu-Mg-Zn* 


Composition 

(Wt%) 

Temperature 

(°C) 

AHf 

(cal/g) 

52% Zn 

340 

43 

— 

588 

55 

56% Si 

946 

181 

20% Si 

803 

47 

16% Ca 

790 

65 

— 

841 

58 

9% P 

715 

32 

30%Cu 5%Si 

571 

101 

25%Cu 15%Zn 

452 

61 

45%Cu 6%Mg 

705 

42 

35%Mg 6%Zn 

443 

74 

28%Ca 17%Zn 

400 

35 

27%Si 17%Mg 

770 

101 

25%Cu 15%Ca 

453 

44 

38%Si 15%Zn 

800 

75 

17%Zn 14%P 

720 

88 

iO%P 7%Si 

840 

22 

19%Zn 7%Si 

765 

30 

30%Mg 21%Ca 

865 

73 

22%Cu 18%Mg 

6%Zn 520 

73 


Density AHf Cp 

(g/cmM (cal/cm^) {cal/g°C) 


— 

— 



1.90 

343 

0.19 

solid 

6.60 

310 

0.12 

solid 

1.38 

90 



5.60 

179 



2.73 

2.80 

276 

171 

r0.31 

^0.29 

solid 

liquid 

8.67 

364 

0.10 

solid 

2.38 

2.26 

176 

79 

rO.39 

^0.35 

solid 

liquid 

4.15 

419 

0.18 

solid 

2.00 

88 



7.00 

616 

0.13 

solid 

6.88 

151 



7.17 

215 



2.25 

164 



3.14 

229 

rO.36 

^0.27 

solid 

liquid 


* Best alloys for heat storage purposes 



the volumetric storage densities were measured with a 
pycnometer bottle. 
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DISCUSSION 


The search method exceeded expectations by 
converging in three steps in the most complicated 
case. SEM metallography was used in almost all 
cases. Even with electron microscopy it was not 
possible to distinguish easily all the phases present. 

In a number of systems the phases had to be recognized 
either by composition or X-ray diffraction measurements. 
Eutectic mixtures also were difficult to recognize; 
lamellar morphologies seldom occurred. In addition to 
metallography, comparison of the overall composition to 
that of the fine, clearly eutectic region was helpful. 
Size of the phases was not by itself an adequate 
criterion. Extreme cases were found. Mg-Cu-Ca had such 
a small eutectic spacing that X-ray diffraction studies 
were necessary to determine whether there was a eutectic 
or a single intermetallic compound. Several cases 
exhibited widely varied spacings in different regions 
of the sample, probably as a consequence of different 
freezing rates. 

Compositions by EDAX can be measured within 3% 
absolute error. The most serious limitation of the 
EDAX measurement is its spatial resolution. If a small 
particle of a phase to be analyzed is embedded in a 
different matrix, the penetration of the beam at the 
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lowest possible voltage in some cases excites the matrix 
elements beneath the particle. Also it is impossible to 
get precise composition values for the different phases 
that constitute an eutectic if the spacing is small. 

The performance of the DTA for large samples 
was entirely adequate for qualitative measurements. 

A 30°C difference between the measured and actual 
transition temperatures was found consistently for all 
the runs owing to the thermocouple position. However, 
the interpretation of the DTA results for multicomponent 
alloys was not always clear. For example, in the 
Al“CU"Si“Mg case four peaks were obtained, one peak 
corresponding to each of the intermetallic compounds 
present in the microstructure. The temperatures of 
these peaks did not correspond to the melting temperatures 
of the compounds that were previously known. 

From an experimental point of view the 
characterization of congruently melting alloys as heat 
storage materials poses no great problems. The commercial 
calorimeters performed well within the expected limits, 
uncertainties of 5% for DTA and 3% for DSC. Conversion 
of heat storage densities from a mass to a volume basis 
requires improved density measurements. The present 
uncertainty in the room-temperatures densities is due mainly 
to the porosity of the alloys. The use of the induction 
furnace, with its vigorous stirring yielded alloys with low 
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porosity. Measurements of densities at the transition 
temperatures are being performed in a related study (20) . 

Although the DCS gave reliable heat capacities for 
alloys that melt below 700“C, for the higher temperature 
alloys the measurement of heat capacities of the liquid 
should be done by some method other than DTA. Vaporization 
losses and reaction of the vapors with the Pt cups are 
likely to be responsible for the failure of DTA to yield 
acceptable heat capacities. The vaporization problem 
does not affect the measurements of heat of fusion 
appreciably because this shift occurs to the baseline as 
a whole. However, there may be a residual effect that 
cannot be determined from the present results. 

As reported previously (19) , the heats of trans- 
formation are systematically well below the calculated 
values as summarized in Table 5. The predictions of the 
ideal model deviate to a large extent from the true 
eutectic temperatures, but the compositions are a useful 
approximation. The regular solution model improves the 
values obtained for the temperatures, but not appreciably 
those for the compositions. In the ternary case there is 
no improvement whatsoever, due to the assumptions that 
had to be made in estimating the ternary activities 
from binary data. 
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Further alloying of the promising binary 
eutectics lowers the transformation temperature as 
expected, but it does not necessarily improve the heat 
storage density as expected. The deviation of the 
actual values from those calculated is erratic and 
unpredictable. The main reason for this variation 
probably is the different degrees of ordering that take 
place in the liquids. Since the presence of stable 
intermetallic compounds was sought to increase the heat of 
transformation, the probability of ordering in the liquid 
is high. To take this complication into account 
additional thermodynamic data are needed. 

The fact that in many cases the addition of a third 
component did not improve the heat of transformation as 
expected suggests that ordering in the liquid is more 
likely for multicomponent systems than for binary systems. 
Statistically it is more likely that some pair of elements 
will exhibit strong affinity in the liquid and therefore 
contribute to short range ordering. 

The practical implication of these conclusions is 
that it will be difficult to find alloys among these metals 
that have storage densities substantially greater than the 
ones presently reported, although some aspects of their 
performance such as volume change during melting may 
continue to improve. Investigation of quaternary alloys 
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may put metals into a more competitive position for 
heat storage at lower temperatures. Further alloying 
may also yield a better accommodation of the melting 
temperatures to the sources. 
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X-ray Diffraction Data for the 
New Intermetallic Compounds 


System 


d (A) in 

decreasing 

order of intensity 

Mg-Cu-Ca 
ternary comp. 

2 . 63 

1.40 

2.24 

1.59 

P-Cu-Zn 
ternary comp. 

1.30 

1 

1 

1.85 

2.09 

2.19 

P-Cu-Si 
ternary comp. 

1 

1.97 

1.93 

1.68 

2.12 1.60 1. 

MgGa2Si2 

2.56 

2.23 

1.41 


Al-Cu-Mg-Si 
quat. comp. 

2 . 59 

3.60 

2.01 

5.55 


1.13 
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Detsrmination of Density Changes During Melting 
by X-Ray Absorption 

C. ERNEST BIRCHENAI.L. ANDREW J. HARRISON, AND SILVIA N. BALART 

Continuous measurement of X-ray transmission through a sample confined in a rigid holder 
permits thermal expansion coefficients of solids, density change during melting, and thermal 
expansion coefficients of liquids to be determined by a single method as temperature is 
increased. Application to melting of pure metals and binary eutectics is illustrated for 
aluminum, the aluminum-silicon eutectic, and the aluminum-copper eutectic. Parameters 
affecting design and operation of the system are discussed. 


Thermal expansion coefficients and volume 
changes during transformation are measured readily by 
dilatometric methods when reactants and products both 
are liquid, both are solid, or when a discontinuous solid 
forms within a liquid phase. However, it is much more 
difficult to measure by dilatometry the volume change 
when a solid completely transforms to a liquid or the 
reverse. Previously, independent measurements of the 
solid and liquid densities at the melting point were 
required. The density of the solid is best determined by 
the extrapolation of thermal expansion measurements. 
The most accurate dilatometric and X-ray diffraction 
measurements on metals have been reported by Sim- 
mons and Balluffi.' Density measurements on many 
liquid metals have been reported by Lucas^ who used an 
Archimedean bob weight method. 

An alternative method is presented in which density 
changes during expansion and transformation are meas- 
ured by the change in transmitted intensity of a 
monochromatic X-ray beam through a path of known 
length. The procedure is straightforward for the melting 
of a pure substance, a congruently melting phase, or a 
eutectic mixture provided that the solid phases in the 
eutectic mixture are uniformly dispersed on a fine scale. 

The method was developed to measure volume 
changes during transformation of some alloys that are 
being considered as possible high temperature heat 
storage materials, to facilitate design of reliable heat 
exchangers. The experimental method and the initial 
measurements on aluminum and two of its eutectics are 
described. 


sity, ja/p = mass absorption coefficient, p = density 
and l - sample thickness or X-ray path length. 

The mass absorption coefficient jx/p is a constant for 
a particular wavelength of X-ray radiation and for a 
particular absorbing element; it is independent of the 
chemical or physical state of the element and orien- 
tation of crystals. Consequently, for an alloy the 
average mass absorption coefficient is given by the mass 
weighted average for all the constituents. If /, is the 
mass fraction of component / in the alloy, 



During congruent melting or eutectic transformation all 
components undergo the same proportionate change in 
spatial density even if some individual solid phases 
increase in density while others decrease, provided that 
no porosity develops. Combining Eqs. [1] and [2], and 
subtracting the logarithm of the intensity transmitted 
through the solid from that transmitted through the 
liquid 

log. (/,//,) = /(P, - [3] 

Equation [3] can be solved for the density change on 
transformation in terms of the change in transmitted 
X-ray intensity. 

To find the volume change on transformation re- 
quires that, in addition to the density difference p, — p, , 
either the density of the solid or of the liquid at the 
melting point be known. 


THEORY 


The absorption of X-rays by any element is given by 



where / == transmitted intensity, /q = incident inten- 
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_ 1/ Pl - 1/p. ^ P._^ P/. (4 

K I/Pa P/ ‘ 

The solid density at the melting point can be obtained 
from the same experiment if the alloy is free of voids, its 
density at room temperature is known, and the sample 
deforms on expansion so that it maintains a thickness 
defined precisely by the rigid cell walls. An analog of 
Eq. [3], including a correction for the thermal expansion 
of the cell from room temperature to the melting point, 
permits evaluation of the linear coefficient k,. of thermal 
expansion of the material. 

= {Prlr - Pm<r [• + k^{T „ ~ 7 ', ) ] } 
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The subscripts r and m denote room temperature and 
melting temperature, respectively. Above the melting 
temperature, the volumetric expansion coefficient of the 
liquid metal can be obtained in a similar way. 

SELECTION OF MATERIALS AND CONDITIONS 

Using the equations presented above, the tabulated 
mass absorption coefficients (Table 1) and densities, the 
feasibility of using this method for the determination of 
volume changes on transformation of the metals and 
metal alloys that are being considered as possible heat 
storage materials was investigated. 

Absorption in the Sample Cell 

The sample cell materials must have a low absorption 
cross section to X-rays (/.e., low atomic number) and 
also have suitable refractory properties. Boron nitride, 
graphite and alumina were considered as candidate 
materials. The fraction of the incident X-ray intensity 
transmitted through the cell walls was calculated using 
Eq. [1]. The fraction of the incident X-rays transmitted 
was calculated for Cu, Mo, Ag and W AT„ radiations 
(Table II) based on a total wall thickness of 4 mm, 
assuming that 2 mm was about the minimum thickness 
per wall for mechanical stability. 


Table 1. 

Mass Absorption Coefficients 
X-Rays (In em’/g) 

* for 

Characteristic 


Cu K, 

Mo K. 

Ag/f. 

w/c„ 

A » 

1.542A 

0.711 

0.561 

0.209 

B 

2.4 

0.39 

0.280 

0.17 

C 

4,6 

0.63 

0.370 

0.18 

N 

7.5 

0.92 

0.520 

0.19 

0 

11.5 

1.4 

0.690 

0.19 

Mg 

39.0 

4.2 

2.117 

0.27 

Al 

48.0 

5.2 

2.64 

0.29 

Si 

62.0 

6.6 

3.26 

(0.33) 

P 

76.0 

8.1 

3.95 

(0.37) 

S 

90.0 

10.0 

4.73 

0.42 

Ca 

65.0 

20.0 

8.82 

0.65 

Cu 

53,0 

51.0 

24.38 

1.8 

Zn 

61.0 

55.0 

26.75 

1.95 

Pb 

250.0 

130.0 

68.08 

5.5 


*These value.s are averaged for internal consistency from tables 
given by C. S. Barrett and T. B. Massalski “Structure of Metals,” 
3rd ed., McGraw-Hill, 1966; L. S. Birks, “Electron Probe 
Microanalysis,” Interscience Publishers, Inc., 1963; and R. D. 
Dewey, R. S. Mapcs and T. Reynolds, “A Study of X-ray Mass 
Absorption Coefficients with Tables of Coefficients,” Reynolds 
Metals Co., Met. Research Div., Richmond, VA (1967), The values 
in parentheses were interpolated. High precision in these coeffi- 
cients is not critical to the design calculations or to the experimental 
results. 


Table II. Fraction of Incident X-Rays Transmitted Through 4 mm 
of Several Low Atomic Number Refractory Materials 



Cu K, 

Mo 

Ag/f. 

W A. 

Graphite 

0.016 

0.567 

0.715 

0.850 

Boron Nitride 

0.012 

0.555 

0.697 

0.850 

Alumina 

v.low 

0.010 

0.112 

0.695 
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This calculation suggests that Cu AT, radiation should 
not be used with any of these refractory materials and 
that Mo should not be used with alumina or a 
refractory material of higher atomic number. With the 
last exception, radiations of Mo, Ag, or W should 
be satisfactory with any of the three materials, including 
thicker cells. 

Absorption in the Sample 

The sample thickness must be large enough that 
errors due to heterogeneities in the microstructure do 
not arise and so that any mismatch between the cell 
walls and the sample is negligible. On this basis the 
minimum sample thickness was taken to be 2 mm. The 
maximum sample thickness is limited by the need for a 
transmitted intensity sufficiently high to give statisti- 
cally reliable count rates. Since the characteristics of the 
X-ray source and of the detector were unknown, this 
thickness was left as an unspecified parameter. How- 
ever, as a first approximation, a transmitted intensity of 
10"^/() was taken to be the minimum for useful 
measurements. This estimate proved to be very con- 
servative. Less than 10~‘/(, was adequate, even when 
the generator was operated at low power. 

Table III lists the fraction of the incident radiation 
that is transmitted through the cell and specimen for six 
combinations of average mass absorption coefficient 
and density to cover the range of likely interest (com- 
pare Table I). The intensity transmitted through the cell 
is taken to be 0.6 /o and the sample thickness 3 mm. 
Cases a, b and c are appropriate for Mo X, radiation; 
cases d, e and f are appropriate for W X, radiation. 

These calculations suggest that Mo X. radiation 
would be suitable only for the light elements, e.g. Mg, 

A1 and Si, whereas W X, radiation could be used with 
any of the elements. Ag X„ radiation would give inter- 
mediate behavior, but closer to Mo X„. 

As Eq. [3] shows, another consideration in choosing 
the appropriate radiation is that the lighter elements 
give a smaller change in signal for a given change in 
density with hard radiation than with soft radiation. In 
order to get the maximum signal change, and with it the 
most accurate determination of the density change, the 
lowest energy radiation that gives sufficient transmitted 
intensity should be used. The last line gives the ratio of 
intensities at the transformation temperature for the 
cases considered if the liquid density is 0.95 times the 


Table III. Fraction ol Transmitted Radiation Through the Cell 
and Sample (3 mm thick) when the Cell Absorbs 0.4 7, 



a 

b 

c 

d 

c 

f 

(^)ave 

5 

15 

30 

0.27 

0.5 

1.5 

\Pf 

P 

3 

3 

5 

3 

2 

5 

//!„ 

6.7 X 10-' 

7.4 X 10--' 

V.low 

0.47 

0.44 

0.063 

/,//, 

1.25 

1.57 

9.5 

1.01 

1.02 

1.12 


/|//, is the ratio of intensities for the transformation temperature if 
the liquid is S pet less dense than the solid. 
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solid density. The small intensity changes for cases d 
and e make it clear that W would be useful only with 
heavier elements or thicker samples. 

EXPERIMIINTAL PROCEDURE 

Description of the Apparatus 

The apparatus is shown schematically in Fig. 1. The 
axis of the windows of the graphite sample cell is 
aligned with a monochromatic X-ray beam. Transmit- 
ted intensity is measured with a Gc detector and an 
energy dispersive pulse counting system. 

The cell. Fig. 2, was clamped firmly between tsvo 
supporting stainless steel plates to compress the sample 
slightly on assembly, ensuring a tight fit. The cell walls 
are thinned in the center to 1.5 mm to minimize 
absorption of the X-ray beam while still retaining 
structural stability. Differently shaped sample retainer 
plates, (cl) and (b) in Fig. 2, had a recess to receive the 
thermocouple. 

The thickness is controlled accurately to better than 
0.2 pet to ensure a tight fit inside the cell. Smaller 
sample sizes could be used if necessary, the main 
limitations on minimum size being the area of the 
sampling beam and the accuracy of the alignment. 
When larger samples of material are available, it is 
practical to use a large beam area to minimize errors 
that could result from small scale heterogeneities in the 
test material and small changes in alignment of the 
system as it heats. The alloy samples are examined 
metallographically to ensure that the phases are dis- 
persed evenly on a scale fine enough for the meas- 
urements to be representative. For example, the Cu-Al 
eutectic alloy measured in this study had an average 
interlamellar spacing of about 3 X 10“’ mm, while the 
beam as it passed through the sample had a diameter of 
about 7 mm. Consequently, the density measured 
should have been a very good average. 


X-RAY 



Fig. I — Schematic diagram for the X-ray absorption apparatus for 
measuring changes in density. 



Fig. 2 — Graphite sample cells with stainless steel supports. 


The vacuum furnace assembly is shown in Fig. 3. It is 
mounted on an aluminum base under a stainless steel 
cover. Beryllium entry and exit windows for the X-ray 
beam 5/8 in. diam and 0.01 in. thick are held between 
0-rings in standard fittings soldered to the cover. A 
rough vacuum, less than 5 torr, was provided to reduce 
absorption by hot gases of varying density along the 
X-ray path. The maximum temperature attainable by 
the furnace is near 1000 °C, but the base and the 
alignment track become excessively hot when the fur- 
nace reaches 800 °C. The fibrous insulation that was 
used initially between the furnace and the stainless steel 
cover gave off a volatile binder that caused tempera- 
ture-dependent absorption by the gaseous impurities. It 
was replaced by two nested stainless steel heat shields. 

From Mo and Ag tubes powered by a high stability 
X-ray high voltage generator. Mo and Ag beams 
are obtained by filtering the radiation with a graphite 
monochromator. The radiation intensity was attenuated 
by Zr foils to avoid overloading the detector. 

Initially, the X-ray output, without passing through a 
monochromator, was metered by a scintillation detector 
in conjunction with a single channel pulse-height ana- 
lyzer. This detector did not satisfactorily resolve the 
radiation and the higher energy continuum radiation 
that is not absorbed very effectively by the zirconium 
filter. The latter tended to overload the detector. To 
overcome this problem, a Ge detector with a multi- 
channel energy-dispersive analyzer provided adequate 
resolution of the K„ peak and integration of this or any 
other selected region of the spectrum. The electronics 
automatically compensated for dead time at high count- 
ing rales. Introduction of the graphite monochromator 
eliminates the need for a beta absorption filter. Also, 
detection could be done adequately with a much 
simpler counter. 

Procedure for Measurement 

A rectangular parallelepiped or circular disk of the 
sample of measured size is weighed to determine the 
sample density at room temperature to about 4 kg/m’. 
The sample, compressed into the cell, is placed inside 
the furnace’. X-ray alignment is verified using a fluo- 
rescent screen. The generator output is adjusted so that 
the count rate is > ICP/s, but so that the ‘dead time’ of 
the analyzer is less than 10 pet. The absorption coef- 
ficient of each alloy is determined by recording the 



Fig. 3 — Vacuum turnacc assembly. 
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reduction in transmitted intensity when samples of the 
alloy of known thickness are used to attenuate the 
beam, Eq. [1], 

The count rate over 100 s intervals is checked for 
repeatability, and the temperature of the sample is 
recorded. The furnace is switched on, and 100 s counts 
are recorded as the sample heats up. The temperature of 
the sample is recorded at the beginning and end of each 
count period. The heating is maintained at about 8 
K/min by manually increasing the furnace supply 
voltage as required to compensate for the increasing 
heat loss with temperature. Near the transformation 
temperature the heating rate is reduced to permit a 
better determination to be made of the intensity change 
on transformation. Measurements are continued above 
the melting temperature to determine the volumetric 
expansion coefficient of the liquid alloy. The maximum 
temperature that can be used is determined by the 
furnace limit (1000 °C), by the volatility of the sample, 
or by misalignment from overheating the base. 

A calibration experiment also is carried out in which 
the cell is empty while it is heated, but the X-rays are 
attenuated by an equivalent sample mounted outside 
the furnace. Thus correction can be made to the 
expansion coefficients for changes in transmitted in- 
tensity that are not caused by density changes in the 
sample. Figure 4 shows that no correction was found 
when using the round cell, but some increase in the 
intensity appeared during heating when using the 
parallelepiped cell. 


RESULTS 

Table IV summarizes the data and results obtained. 
The plot of logarithm of transmitted intensity V5 average 
temperature for A1 is shown in Fig. 5, for typical 
experiments. The expansion data for the solid and for 
the liquid were fitted with a third order polynomial 
expression using a conveniently available computer 
program. The values of the transmitted intensity at the 
transformation temperature for the solid and the liquid 
were determined from extrapolation of these fitted 
curves. 

The change in transmitted intensity at the trans- 
formation temperature was reproducible to about 0.1 
pet. However, for solid expansion near room temper- 
ature the apparent coefficient of thermal expansion was 
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Fig. 4 — Calibration of empty cells show that no correction is required 
for the round sample holder. Correction is necessary above 500 ‘^C for 
the parallelepiped cell. 


not always reproducible. This difficulty is not com- 
pletely resolved. If some specimens have undetected 
porosity at room temperature, they should exhibit 
abnormally low coefficients of solid expansion over 
some range if the pores collapse during heating below 
the transformation temperature. If the pores persist to 
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Fig. 5 — Plot of log intensity against temperature for commercially 
pure (1 100 Alloy) and research grade (99.99+ pet) A1 illustrating 
typical data. The broad wave above 400 *C, typical of the com- 
mercially pure A1 but not of high purity Al, is not reproducible in 
detail. The intensity axis displays only relative values. An arbitrary 
shift was used to avoid overlapping points. The samples did not have 
the same thickness. 


Table IV. Density Changes for Al and Two Al-Eutectic Alloys 


Density 

in kg/m’ 

AK/F,t 

T'C 

Measured 

Tabulated* Measured Tabulated* 



Pure Al 


20(s) 

2690 

2699 (6) 

0.0769 (3) 

660(r) 

2558 

2555(1) 

0.072 0.0759 (4) 

660(f) 

2377 

2368-2410(6) 


760(f) 

2343 

2338 (6) 

0.069 (5) 



Al-17.5 at. pet Cu 


20(.r) 

3506 

3570 ) 


548 (s) 

3424 

3400 1,,, 

0.051 0.050 (5) 

548(f) 

3258 

3290 

748(f) 

3186 

3200 ) 




Al-12 at. pet Si 


20(r) 

2686 

2693 (5) 


579 (i) 

2553 


0.048 0.040 (5) 

579(f) 

2445 



779(f) 

2382 



•Numbers in parentheses cite references. 
tAt the transformation temperature. 
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the transformation temperature, the solid density re- 
ported at that temperature would be correspondingly 
low . The expansion coefficients reported in Table V are 
from those runs that had no unusual trends in the solid 
Expansion curve. 

ii Pure aluminum. The volume change on melting was 
found to be 0.072 of the solid volume at the melting 
point. This value is an average from those experiments 
which gave reasonably good solid expansion data. 

Single results were between 0.071 and 0.073. The figure 
is a little lower than the value 0.0769 determined by 
3\'ray’ from estimates of the best density measurements 
on the solid' and liquid.^ However, they are greater than 
■ the value 0.069 reported by other workers before 1970.* 

ii) Al-12 at. pet Si. This eutectic alloy showed a lower 
change in volume than pure aluminum. The value 0.048 
gisen in Table IV is the average for two experiments. 
Although this is a commercially important system, 
precise values could not be found in the literature. The 
tabulated value reported here is interpolated between 
two values reported by Mondolfo.* 

Hi) A 1-1 7. 5 at. pet Cu. This alloy, measured before 
installing the monochromator, is unfavorable for Mo 
radiation because of the substantial Cu content, which 
strongly absorbs the radiation. A narrow band of the 
higher energy white radiation was selected by means of 
the energy dispersive analyzer, to monitor the density 
changes. The result, a change of 0.051 of the solid 
volume, is in excellent agreement with the literature. 

For the three kinds of samples, the expansion data 
exhibit greater scatter for the liquids than for the solids 
because only about 10 data points are obtained before 
the cell base overheats. This shortage of data con- 
tributes most of the error in the final result. 

DISCUSSION 

This X-ray absorption method should be applicable 
to the measurement of coefficients of expansion and 
volume change during transformation of any material 
that can be prepared in non-porous condition and that 
also has sufficient ductility to flow when heated inside 
the absorption cell. To avoid porosity after casting, the 


Table V. Linear (Solid) and Volumetric (Liquid) Expansion 
Coefficients 

Linear Expansion Coeff. Volumetric Expansion Coeff. 

Measured Tabulated* Measured Tabulated* 

Pure Al 

2.77 X 10-’ 2.95 X I0->(2) 0.9 x 10-* l.l X 10"* (6) 

2.50 X 10- ’(6) 

range: 20 to 660 *C range; 660 to 760 *C 

Al-12 at. pet Si 

1.7 X 10-* 1.3 X I0-* 

range; 20 to 579 *C range: 579 to 679 *C 


Al-17.5 at. pet Cu 

1.5 X 10 ’ 1.1 X 10-* 

range: 20 to 548 *C range; 548 to 748 *C 

* Numbers in parentheses are references. 
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alloys were rcmelted at moderate pressure,’ yielding 
good results. However, it is possible that some of the 
discrepancies remaining in the lower temperature ex- 
pansion results are due to squeezing out undetected 
microporosity as the solid is heated in the absorption 
cell. 

In the case of pure elements the undesirable effect of 
impurities is clearly shown during expansion. For 
commercially pure Al (1 100 Alloy), a bump of variable 
size was found in the expansion curve between 400 and 
500 °C which did not appear when using research purity 
(99.99 pet) Alcoa aluminum. However, the two grades 
had indistinguishable volume change during melting. 

When a density change is determined by taking the 
differenee between the extrapolated density of solid and 
liquid at the transformation temperature, the uncer- 
tainties in the solid and liquid densities add to become 
the uncertainty in the small difference. Very careful 
measurements are required on both states of each 
system. Perhaps this is the reason that reliable volume 
changes are available for melting of only some of the 
elements and few of the binary or higher eutectics. 
Nevertheless, some regularities are discernible among 
these results, probably because the liquids share a 
nearly common type of packing regardless of the stable 
crystalline structures that melt. 

Tabulated values' ** show that the cubic close-packed 
metals usually contract about 5 to 5.5 pet relative to the 
liquid volume when they freeze. Lead appears to 
contract a little less (3.6 pet) and aluminum a little more 
(7.2 pet), as confirmed for the latter in this study. The 
hexagonal close-packed metals cadmium and zinc con- 
tract about 4 pet upon freezing. The less dense solid 
structures, such as bcc thallium and iron, and the 
tetragonal elements tin and indium contract only about 
2 to 3 pet. The rhombohedral, layered semimetals 
bismuth and antimony, and the covalent semiconduc- 
tors silicon and germanium, that crystallize in the open 
diamond cubic structure, expand when they freeze. 

A eutectic mixture should exhibit a volume change 
upon congruent melting that is the volume-weighted 
sum for the constituent solid phases. Therefore the 
aluminum-silicon eutectic should have a smaller in- 
crease in volume than does pure aluminum because the 
silicon probably has a smaller partial molar volume in 
the melt than it has as a nearly pure constituent in the 
solid eutectic mixture. The similarly smaller volume 
change for the aluminum-copper eutectic may be a 
consequence of the relatively open structure of the 
intermediate 0 phase of composition AIjCu. 

CONCLUSIONS 

The X-ray absorption method appears to be capable 
of determining density changes during melting with 
higher precision than taking the density differences 
from two different types of measurements, one on the 
solid and another on the liquid. Small volume changes 
might be achieved if intermediate phases with open 
structures that result from strongly covalent or ionic 
contributions to binding are part of the mixture in 
sufficient concentration to counteract the expansion 
upon melting of the metallic solid solutions. 
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APPENDIX D 


VOLUME CHANGE DURING MELTING 
OF EUTECTIC ALLOYS 


BRUCE PREGGER 




Radiation Source 

In previous work the source of radiation used 
was an X-ray tube powered by a Philips XRG-3000 high 
voltage generator (Appendix C) . An Ag tube was used to 
provide AgKa X-rays at 22Kev. The high intensity of this 
source necessitated the use of a graphite monochromator 
and Zr foil filters. Results for the volume expansion of 
Al, Al-0.12Si, and Al-0.18Cu were obtained with this 
system that agreed with previously reported data. However, 
several characteristics of this source were undesirable. 
Variations in the observed intensity were systematic in 
nature and were attributed to a real variation in the 
output of the X-ray tube over the course of an 
experimental run. These variations were dependent on the 
change of the rate of flow of cooling water to the generator, 
an^ possibly also on variations in line voltage. The total 
intensity of the tube was a problem. The simple addition 
of absorbing Zr foils was not satisfactory, because high 
energy continuum radiation was relatively unaffected and 
remained to saturate the detector. A monochromator 
eliminated the unwanted radiation, though problems in beam 
alignment resulted. After alignment, however, the 
monochromator-Zi* foil combination produced an intensity 
which yielded good counting statistics with an intrinsic Ge 
detector. Finally, the nature of the experiment required 
that the sample furnace be placed in close contact with 


Dl 



the X-ray generator, a situation which led to the premature 
termination of experimental runs as the furnace became hot. 

A theoretical limit of 1000°C was obtainable with the 
furnace, but at about 800°C the metal base reached what 
was determined to be a maximum safe temperature. 

The use of a radioactive isotope as an X-ray source 
in place of the X-ray tube and generator offers a number 
of advantages. The output intensity of an isotope is 
constant over experimental times within the limits of 
radioactive decay statistics. Such a source is small and 
easily managed; alignment can be facilitated. The source 
can be physically separated from the furnace, allowing 
maximum temperatures. The radiation from an isotope is 
characterized by sharp energy peaks, so that the use of an 
energy dispersive counting system with such a source 
makes possible measurements without a monochromator. The 
judicious choice of an isotope, most of which have more 
than one useful energy peak, allows measurements to be made 
on materials of very different densities. From these 
considerations it was concluded that the volume change studies 
should continue with an isotope replacing the X-ray tube/ 
generator . 

The choice of an isotope was dependent upon the 
following factors — the energy peak (or peaks) had to 
optimize the volume change measurements for the greatest 
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number of alloys, and the isotope had to be commercially 
available at a reasonable price. 

To optimize volume change measurements the 
equation 

2(/N +/BKG) 

N = number of counts - background (BKG) 

2a = uncertainty/measurement 

must be maximized. This depends upon the thickness of the 
sample, a controllable parameter, the number of counts per 
measurement and background, which are somewhat adjustable, 
and on the volume change and mass absorption coefficient, 
determined by the material. 

The alloys to be measured fell into two groups 
separated by their absorption characteristics — those 
containing Cu or Zn, and those that did not. Cu and Zn 
absorb radiation much more strongly than Al, Si, Mg, or 
P, and thus an alloy with either one must be measured with 
high energy X-rays. Radiation of this energy does not 
show enough variation with density to be of use with the 
lighter alloys. Another consideration is the detector, 
as a Nal crystal and a Ge detector behave considerably 
differently at different energies. 
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Calculations were done to compare the relative value 


of two promising isotopes, Cd-109 and Gd-153. Cd-109 

has peaks at 22, 25 and 88 KeV, Gd-153 at 41 and 47 KeV. 

Calculations were made for the 22 KeV Cd and the 41 KeV 

Gd peak as they would be detected with a Ge solid state 

detector, and for the 88KeV Cd peak as it would be 

detected in a Nal crystal. A series of calculations was 

performed for each alloy at each wavelength in order to 

determine the sample thickness that would give the best 

volume change precision. Values of density were measured 

from our prepared samples. Mass absorption coefficients 

were estimated by summing the mass fractions of the 

* 

coefficients of the component elements. Elemental 
absorption coefficients are given in Table 1, estimated 
alloy coefficients in Table 2, and optimum volume change 
precisions with corresponding sample thicknesses in Table 3. 


From Table 3, it can be i^een that the 41 KeV Gd peak 
with Ge detector gave the best precision for the dense 
Cu,Zn alloys. For the lighter alloys the 22 KeV Cd peak 
gave the best results. This peak energy could not be used 


* 

Volume change estimates were made by summing the mass 
fractions of reported elemental change values. The count 
rate was estimated through consideration of detector 
efficiencies and collimated radiation beam intensities 
produced by 5 mCi sources. Counting time was taken to be 
1000 seconds. 
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Table 1 


Mass Absorption Coefficients for Cd-109 
Gd-153 (Cm^/gm) 


KeV 

Mg 

A1 

Si 

P 

Cu 

Zn 

In 

22 (Cd) 

2.12 

2.60 

2.7 

3.95 

24.38 

26.75 


41 (Gd) 

0.4* 

0.5* 

0.6* 

0.6* 

4.0* 

5.0* 


88 (Cd) 

0.17 

0.19 

0.2* 

0.2* 

0.67 

0.78 

2.2 


Sources: "Handbook of Chemistry and Physics," 50th edition, 

C. R. C. Press, 1970; C. S. Barrett and T. B. Massalski 
"Structure of Metals," 3rd edition, McGraw-Hill, 


1966. Starred values are interpolated. 



Table 2 


o 


Estimated Mass Absorption Coefficients for Alloys 


KeV 

Al-0.12Si 

Al-0.12Si-0.05Mg 

Si-0.35Cu-0.28Mg 

Cu-0.25P-0.14Zn 

Zn-0. 40CU-0. 

22 (Cd) 

2.6 

2.5 

10.3 

20.6 

24.3 

41 (Gd) 

0.5 

0.5 

1.0 

3.0 

3.0 

88 (Cd) 

0.19 

0.19 

0.36 

0.57 

0.65 



Table 3 


Calculated Relative Precisions for Volume Change 
Measurements; (optimum sample thicknesses (cm)) 

KeV Al-0.12Si Al-0.12Si-0.05Mg Si-0 . 35Cu-0 . 28Mg Cu-0.25P-0.14Zn Zn-0.40Cu-0.15Mg 

22(Cd) 8. 5(0. 4) 8. 1(0. 4) 

41(Gd) 4. 8(0. 5) 4. 6(0. 5) 1.1(0. 3) 9. 3(0. 2) 6. 8(0. 2) 


88 (Cd) 


0.26(0.5) 


3.4(0. 5) 


2 . 6 ( 0 . 5 ) 



with the Cu and Zn alloys, but the 88KeV peak can, 
though with less precision than the Gd peak. It was 
decided to use the Cd-109 source. The better results 
for the light alloys would be useful for comparison to 
the results obtained for such materials with the X-ray 
tube used previously. 

APPARATUS 

Radiation Source 

The Cd-109 isotope used was in the form of a small 
disk, 0.79 cm in diameter and 0.51 cm in thickness. It 
was supplied by New England Nuclear (model NER-465 LE66) 
and had an activity of 5 mCi. The source is collimated 
into a narrow beam by having it rest at the closed end 
of a 15 cm stainless steel tube, held in place by a 
smaller diameter tube of the same material which slides 
into the first tube. The second tube has an inside 
diameter of 0.4 cm which acts as the collimator. The 
0.4 cm diameter corresponds to the active area of the disk 
source, ensuring the greatest possible radiation intensity 
for the collimated beam. These tubes containing the 
source are sunk into a lead brick, as shown in Figure 1. 

Sample Furnace 

The furnace is basically that used by Harrison 
(Appendix C) with some minor modifications. Five meters of 
Nichrome wire are wound about the outside of a 9 cm diameter 
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X-RAY SOURCE holder 
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alumina tube. This tube has two 1 cm diameter holes 
drilled through it in order to let the X-ray beam pass. 

The tube sits upright in the center of the furnace on a 
refractory base. The sample holder is placed inside 
this tube. The furnace is designed to operate under 
vacuum. This vacuum is maintained in a stainless steel 
cover which fits onto an 0-ring in the furnace base, and 
which has two windows built into it to allow the entrance 
and exit of the X-ray beam. These windows are made of 
beryllium, 0.01" thick, held in place by 0-rings in 
stainless steel fittings soldered to the cover with silver 
solder. Thermocouple wires, electrical connections, and 
the vacuum line penetrate through the base of the furnace. 
The thermocouple is Chromel-Alumel . A vacuum of less than 
5 torr is realized with a small mechanical vacuum pump. 
Between the furnace windings and the vacuum cover are 
two stainless steel heat shields and 1 cm thickness of 
heavy Alundum refractory. The furnace is shown in Figure 2. 

Radiation Detectors 

Two types of radiation detectors were used, depending 
on the energy of the radiation being measured. For low 
energies, the Cd-109 22 and 25 KeV peaks, a General 
Electric 400 high purity germanium spectrometer was found 
to give the best results. The resolution of this device 
at 22 KeV was measured to be 2.3%. At this energy a Nal 
scintillation crystal had very poor resolution, and could 
not resolve the 22 and 25 KeV peaks and thus could not 
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measure truly monochromatic radiation, leading to 
erroneous results. However, as the energy increases the 
efficiency of the Ge counter drops rapidly. At 88 keV 
the counting rate of the Ge detector is at such a level 
that a Nal crystal is of more value as a monitoring 
device, even though the background count rate in the 
scintillation crystal is much greater. 

The apparatus is arranged on a nonconducting 
bench top in such a manner as to allow changing from 
the solid state to the scintillation counter. The isotope 
holder and furnace are both small and easily repositioned, 
being physically separate. 

Either detector was connected to a Tennelec TC 203BLR 
linear amplifier, from which the signal was run to a 
General Electric 700 multi-channel analyzer. 

The arrangement of the apparatus is given in 
Figure 3. 

Sample Holder 

The sample holder has three requirements. It must 
keep the sample thickness constant throughout the 
experimental temperature range; it must allow enough 
radiation through the sample to provide good counting 
statistics; it must be made of material that will not 
react with the samples at elevated temperature. Graphite 
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was chosen as the sample containment material. Three 
graphite plates are placed together, the center plate 
cut out in the center to contain the sample. The outer 
plates are thinned to 0.2 cm in their centers to allow 
passage of the radiation. These three plates are 
sandwiched between stainless steel plates and pressed 
together tightly with eight screws. The steel plates 
have holes cut out of their centers to allow X-ray 
transmission. The configuration is shown in Figure 4, 
Graphite was chosen as containment material because of 
ease of fabrication and low X-ray absorption character- 
istics. The 0.4 cm total path length for X-ray absorption 
in the sample holder is calculated to transmit 90.9% of 
incident radiation at the 22KeV Cd-109 peak, and 95.8% 
at the 88 KeV peak. 

To facilitate alignment of the sample holddr with the 
furnace windows, radiation source, and radiation detector, 
the holes in the stainless steel plates were threaded 
and fitted with 5 cm long tubes of the same material. 

These tubes defined the radiation path through the 
furnace winding and sample holder. 

Samples 

The aluminum and Al-0.12 Si samples were commercial 
materials. All other alloys were prepared by heating 
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constituent elements and various commercial Mg, Si, 

P-metal compounds in an Ajax-Northrup induction furnace. 

The preparations were sealed in alumina crucibles to 
prevent vaporization of the volatile elements, and heated 
in an argon atmosphere to prevent oxidation. The induction 
heating provided good mixing of the alloys so that homo- 
geneous eutectic structures resulted. The alloys were 
machined into disks that would fit tightly between the 
two outer graphite plates in the sample holder, but smaller 
in diameter than the hole in the center plate to allow 
for transverse expansion on heating and melting. The 
sample thicknesses were determined by optimum volume change 
calculations (see Table 3) . 

The alloys themselves were those which showed 
promise as heat storage materials. Eutectic compositions 
in the various systems were found by Parkas, Appendix B. 

The heats of fusion were measured in the same work. 

Those alloys which demonstrated the greatest heats of 
fusion have possible commercial utility and were chosen 
for this study. Their compositions are given in Tables 
2 and 3 . 

Experimental Procedure 

The alloy sample is placed inside the sample holder 
which is then screwed together tightly. The sample holder 
is positioned in the furnace so that the radiation passes 
from the source to the detector by way of the radiation 
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path tubes on the sample holder. The radiation from a 
5 mCi Cd-109 source is too weak to permit the use of a 
fluorescent screen to check alignment; accuracy in this 
regard must be determined by counting with the detector. 
After satisfactory positioning of the sample holder in 
the furnace, the furnace cover is put in place and the 
vacuum pump is turned on. With the vacuum at 5 torr or 
less, counting may begin. The integral number of counts in 
the proper peak is recorded by the multi-channel analyzer 
for a period of 1000 seconds. A variable output trans- 
former which controls furnace heating is turned on and 
the furnace is allowed to heat so that the temperature, 
as recorded by the Chromel-Alumel thermocouple, is raised 
by about 100 ®C. At this point the transformer is adjusted 
until a temperature equilibrium is achieved. The multi- 
channel analyzer is turned on to record for another 1000 
seconds. Changes is furnace temperature in this time 
period can be minimized with a skilled hand at the 
transformer; such variations are accounted for by averaging. 
The process is then continued, raising the temperature in 
100° intervals through the melting temperature and into 
the liquid range as far as the furnace allows. Measurements 
are also taken for 1000 seconds during the heating intervals 
for comparative purposes, in which temperatures are recorded 
at the beginning and end of the 1000 sec interval and 
averaged . 
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RESULTS 


The values obtained for volume changes of the 
various alloys studied are given in Table 4 . Figures 5 
and 6 show plots of An intensity vs. temperature for the 
various materials studied. The solid and liquid phase 
data points were separately characterized as straight lines 
by the method of least squares. It was felt that the 

scatter of the data precluded the use of higher order 
polynomials. These lines were extrapolated to the melting 
point of the material. The values of the solid and liquid 
line at the melting point were taken as the natural log 
of the incident intensity produced by the solid and liquid 
melting point densities, and as such were used in calculating 
those densities. The value of the solid phase data line 
was used at room temperature to represent the room 
temperature density; the measured room temperature 
intensity was not. Precision was gained in this manner. 

Numbers obtained for the solid and liquid 
expansion coefficients were not in agreement with 
previously measured values, nor were they to any great 
extent reproducible. This is a serious problem, for it 
can affect the fusion volume change calculation. The 
fusion volume changes calculated appear to agree with those 
that have been determined by Harrison and others (2-4) . 

This study is concerned primarily with new alloys on which 
no data of this nature are available for verification. 
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Table 4 


Volume Changes on Melting of In, Al, and 
Various Eutectic Alloys 

(AV/Vg) Measured Tabulated 


In 

0.031 

0.025 

(1) 

Al 

0.09 

0.072 

(2) 



0.0769 

(3) 


Al-0. 12Si 

0.041 

0.040 

(2) 



0.038 

(4) 


Al-0.12Si-0.05Mg 0.066 

Cu-0.25P-0.14Zn 0.033 

(1) Smithells, "Handbook of Metals" 1955. 

(2) Harrison, Appendix C. 

(3) P. J. Wray, Met. Trans., 1974, vol. 5, p. 2602. 

(4) L. Mondolfo, "Aluminum Alloys," Butterworths , London, 1976 
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Al-0.12Si had a change of volume on melting 
that was measured to be 4.1%, very close to the change 
of 4.0% given by Mondolfo (4) and 4.8% as measured by 
Harrison (2) with this method using an X-ray tube source. 

In, measured with the 88 KeV peak, gave a value 
of 3.1% for the volume change. Smithells (1) gives 
2.5%. 

The runs made on pure aluminum were not altogether 
satisfactory. High volume changes were calculated from them, around 
20% higher than the grouping of other experimental values 
(2-4) . The reason for this remains unknown. 

Al-0.12Si-0.05Mg gave a measured volume change of 
6.6%. This value is probably too high, although there are 
no other measured values to check against. The presence 
of Mg has been found to be detrimental to the acquisition 
of good data when present in larger amounts in a graphite 
container; perhaps it is not helpful in this case. 

Cu-0.25P-0.14Zn is an easy-to-handle alloy. The 
volume change measured was 3.3%, which in concert with a high 
heat storage capacity and ease of handling (it was made 
successfully in a resistance furnace) makes this material 
look very promising for commercial utilization. 
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DISCUSSION 


The variability of solid and liquid expansion 
coefficients demands resolution. The problem stems from 
the changing of the radiation intensity in other parts of 
the furnace than the sample. Due to the low level of 
intensity produced by the radiation source, as wide 
a beam as possible is used. This leads to blocking by 
opaque surfaces in the furnace. When the furnace heats, 
these surfaces may moye. in such a manner as to increase or 
decrease the theoretical incident intensity so that the 
equation for p^p - Pj^QQjyj ^ does not hold. To decrease 
the beam diameter to that usable with an X-ray tube 
source would decrease the intensity from the radioactive 
source to a level that would eliminate any chance of 
useful results. Experiments were done using an empty 
sample holder to test this theory of changing incident 
intensity, and what was found was that the furnace did 
change the intensity on heating, although not always by 
the same amount nor always even in the same direction. 

The Pjj - Pg is relatively unaffected by this problem coming 
from data of a single temperature, but if p is calculated 
wrong then the percentage volume change on melting will 
not be accurate. 

A solution to this problem involves the simultaneous 
measurement of two peaks of different energies. Opaque 
blocking accounting for a changing incident intensity 
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may be eliminated between equations for the absorption 
at the two energies giving the following equation: 


I ** “ I ** 0X0 — 

adjusted measured ^ 


fn"'''’!!. r^2 


JjTFT if I - 




E - (SD 


where 


^"adjusted ~ intensity that would be measured had 

not the incident intensity been changed by 
opaque blocking from T' to T" 
y/Pl 2 = mass absorption coefficients for radiation 
of energies 1 and 2 

= measured intensity at T' with energy 1 (for 
example) 

^"measured ^ ^1 ^2 ' yielding corrections at both 

energies. 


This equation accounts for partial opaque blockage, but 
there is no solution for the opposite effect, partial 
free transmission; Thus any porosity in the sample 
cannot be mathematically eliminated. Samples prepared 
in the induction furnace have shown good density character- 
istics, without excessive porosity. 

Cd-109 has two peaks in the same energy range: 

22 and 25 KeV. Test measurements were made on Al for 
corrections by the above equation. The 25 KeV peak is 
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relatively weak, yielding poor counting statistics, and 
results were inconclusive. Gd-153 on the other hand has 
two strong peaks in the same range, at 41 and 47 KeV. 

An active Gd source could be a solution to the problem. 

Another problem was encountered when dealing 
with alloys rich in Mg. At elevated temperatures the 
Mg reacted with the graphite of the sample holder 
changing both sample thickness and sample composition. 

Such was the case for Si-0 . 3SCu-0 . 28Mg and Zn-0.40Cu-0.15Mg. 
A new material must be found to contain these alloys. 

It appears as though a Gd-153 source could provide 
better results than the present Cd-109. The source 
and experimental configuration now being used do not 
seem wholly satisfactory, especially for expansion 
coefficient determination. 


D25 




APPENDIX E 




ANALYSIS OF CYLINDRICAL HEAT EXCHANGER WITH PHASE CHANGE MATERIALS 


Michael B. Labdon Selguk 1. Giigeri 

Graduate Student Assistant Professor 


Mechanical and Aerospace Engineering Department 

University of Delaware 

Newark, Delaware 19711, U.S.A. 


ABSTRACT 

Within the field of energy, there is an increasing emphasis 
on the use of phase-change materials (PCM) for thermal energy 
storage, but the analytical and numerical methods available to 
describe their thermal behavior in cylindrical coordinates are 
relatively few in number. In addition, these methods are either 
mathematically too complex for engineering purposes, or simplified 
to the point of solving a pseudophase-change process far removed 
from the actual system itself. In this paper, an efficient, nu- 
merical method is presented for modeling one-dimensional heat con- 
duction (including sensible and latent heat) around an internally- 
convected cylinder of PCM. Finite-difference equations monitor 
the heat transfer within the solid and liquid regions, while Tay- 
lor-series expansions of temperature are made in the regions close 
to the interface. Inputs to the routine include temperature- 
dependent properties, initial temperature distributions, time- 
varying convection boundary conditions and time-varying inlet fluid 
temperatures. Outputs from the method include melt-line locations 
and temperature distributions vs. time, as well as charging rates 
and pipe surface temperatures. 


1. INTRODUCTION 

Many heat transfer processes include a change of phase of a 
material due to melting or solidification. Ice formation, abla- 
tion of missile skins and solidification of metal castings are a 
few such applications. Within the fields of Thermal and Solar 
Energy Storage, there is an increasing emphasis on the use of such 
phase-change materials (PCM) , for the collection, storage and trans- 
mission of the forms of energy which today are more costly to gen- 
erate and more actively conserved. The many applications of, and 
increased interest in, PCM’s require accurate mathematical model- 
ing of the thermal behavior of such materials. Predictions of 
melt-line locations and temperature histories throughout the medium 
are useful in: a) choosing the appropriate material for PCM con- 

tainment, b) optimum design of the PCM apparatus, and c) the con- 
trasting of several PCM's to choose those applicable for a specific 
use. Due to the inherent temperature non-linearity at the inter- 
face, mathematical modeling of the phase-change process via an 
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analytic approach is quite complicated. As a result, solutions 
exist for systems with simple geometries and very restrictive 
boundary and initial conditions. The extension to systems in 
cylindrical coordinates further reduces the number of available 
solutions due to the added mathematical complexity of the 1/r 
term in the original heat conduction equation. However, this is 
the system geometry often used with PCM's and energy storage or 
transmission; such as that found with heat exchangers, heat pipes, 
tube flows, etc. 

Authors such as Meyer [2], Foots [3], Solomon [4], Boley [5], 
Krieth [6J, JiJi 17], Shubin [8] and Sikarskie [14] have attempted 
such modeling, but required such restrictions as constant wall 
temperature, constant material properties or temperature distri- 
bution assumptions. The intricate mathematics or system-simplifi- 
cation of such methods preclude their use for engineering purposes. 
This is also evident in the numerical solutions alternately em- 
ployed to model the phase-change process. Talwar [9] and others 
reduced the problem to one of pseudomelting Whereby the thermal 
properties at the interface were "smoothed out" to provide a con- 
tinuous relationship between the solid and liquid regions. Forster 
[10] considered the use of a quadratic temperature profile through 
the melt-line between two known points in an effort to compute tVie 
temperature gradients at the melt-front and hence its velocity. 
Murray and Landis [11] further refined Forsters model by institu- 
ting temperature profiles three nodes from the interface, but re- 
sulted in starting and terminating errors that propagated through- 
out the solution. Lazaridis [12] attempted to expand Murray and 
Landis' work to multidimensional systems, and resulted in develop- 
ing a routine both mathematically and computationally difficult to 
implement with various geometries. 

This paper outlines an efficient, numerical routine for solv- 
ing one-dimensional heat conduction (including sensible and latent 
heat) around an internally-convected cylinder of PCM, by noting 
the dominant mode of heat conduction to be in the radial direction. 
Finite-difference equations monitor the heat conduction in the re- 
gions away from the fusion front, while Taylor-series expansions 
of temperature are made near the melt-line, in such a manner as to 
continually trace its propagation throughout each node, while ad- 
justing all nodal temperatures. The versatility of the method is 
reflected in the fact that temperature-dependent properties and 
initial temperature distributions may be included, as well as time- 
varying boundary conditions. Time-dependent inlet fluid tempera- 
tures can also be handled which represents the more realistic 
condition encountered with energy storage and transfer systems. 


2 . PROBLEM STATEMENT 

The heat exchanger geometry under consideration is depicted 
in Figures 1 and 2. Figure 1 is an axial view of a cylinder sec- 
tion of length "S,", containing a PCM partitioned into "N" equal- 
thickness, concentric, annular regions or nodes, with a represen- 
tative melt-line location. Heat convection occurs at the inner 
radius Ri, and, for the analysis, one-dimensional solidification 
occurs radially outward. The outer radius R 2 f is an insulated 
boundary, but may also be subjected to convection if required. As 
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the convecting fluid tra- 
vels through the cylinder 
section an energy balance, 
described in Section 4.5, 
is performed on its bulk 
temperature to arrive at an 
increased outlet tempera- 
ture. This in turn becomes 
the inlet flow temperature 
for the next downstream 
section. This process is 
repeated for "M" sections 
until the desired pipe 
length of "MJl" is attained, 
as depicted in Figure 2. 

In this manner, by taking 
a sufficiently large number^ 
of axial segments, the heat 
exchanger can be effective- 
ly modeled as a two-dimen- 
sional system with the one- 
dimensional radial heat 
conduction mode dominating 
the axial conduction mode. 



Figure 1. Axial view of heat-exchanger 

with representative concentric 

3. PROBLEM SOLUTION rings (nodes) of PCM, depicting 

a particular melt-line location. 

The only requirements 
of the following analysis 

are that the fusion temperature is constant, the density is con- 
stant in both phases, and that natural convection effects are neg- 
ligible. Use of the foregoing restrictions is acceptable for many 
engineering purposes, while exclusion of the same would result in 
a more complicated system modeling. As previously mentioned, only 
radial heat conduction will be considered. 


The problem analysis is facilitated by the introduction at 


PCM 




SECTION M 


Figure 2. Horizontal view of heat exchanger of length "Mil" contain- 
ing a PCM between an internally-convected pipe and an 
insulated outer one. 
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this point, of the governing equations used throughout the phase- 
change medium: 


1 3T^ , 1 3T^ 

3r^ r 3r 3t 

3^T2 1 3T2 1 3T2 

3r^ r 3r tt 2 3t 

3T, 3T, dS 

k, — - - k, — = LP — 

3r 3r dt 


Rl 1 r £ S 


S £ r £ R2 


r = S(t) 


(la) 


(lb) 


( 2 ) 


T 


1 



r = S(t) 


(3) 


where subscripts 1 and 2 denote the solid and liquid regions re- 
spectively. Eqs. la, b are the field equations in the solid and 
liquid regions. Eq. 2 is the well-known Stefan Condition result- 
ing from an energy balance at the interface, and Eq. 3 is the con- 
dition of constant fusion temperature. As an added convenience, 
the following non-dimensionalizing parameters are introduced into 
the above equations: 


C = r/R ; T t/R^ 


g(t) = S(t)/R ; A? = Ar/R ; K = k 2 /kj^ 


~ “ 2 ^“! ' ^ ~ 


6? s 6r/R ; = R^^ ; = R 2 


where R is a chosen reference length. The resulting equations, in 
which all variables except temperature are non-dimensionalized, ap 
pear as: 


3^Tj^ 1 3Tj^ 3T^ 

3C"' C 3t 



1 ^ 
K 9^ 


9t 


^1 - ^ ^ 
^ 1 ^ 1 


3T, 3T, de 

— - - K — - = X — ? = e(T) 

3^ 3£ dx 


(4a) 


(4b) 


(5) 


Figure 3 illustrates the cross-sectional, radial view of a 


E4 


T 



Figure 3. Radial Temperature distribution near the interface, with- 
in a particular pipe section. 


representative pipe section ”i”, 1 £ i £ M, undergoing solidifica- 
tion at some intermediate time "t", where the melt-line has just 
traveled to the distance ei(T). The node in the solid region 
closest to the interface is denoted as "q", while "q + 1" repre- 
sents the similar node in the liquid region. The nodal spacing, 
ACf is defined as (Cn “ ?l)/(N - 1) while 6C is the time-varying 
distance between node q and ei. At this particular time t, only 
the current location is known, as well as the "old" (t - At) 

nodal temperatures i = 1, M; j = 1, N. The calculation of 

"current" temperatures and future melt-line positions requires the 
successive execution of the following computational logic: 

a) Determine all nodal temperatures except T^^q and 
at time t via the finite-difference equations, 

b) Compute the current temperature gradients at the 5 nter- 

face, 9Tj^/8^| = A^ and 9T2/9^ 1 = and, 

^i ^i 

c) Determine the current values for temperatures T^^q and 
Ti,q+i, as well as the melt-front velocity dci/dx/ 

all of which are outlined in the next three sub-sections. 


3.1. Finite-Difference Equations 

Temperatures at all nodes i,j where 1 ^ i ^ M and 1 £ j _< q-1. 
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q + 2 ■< j •< N are computed from the usual finite-difference approx- 
imation to the original heat conduction equation (4) , 



where Tj,^j denotes the new value for j after the time increment 
At, and y"^= 1 when j is within -{ihe sol i<5"^ region or y = Yjj, when j 
is within the liquid region. is solved from (6) as. 


1 

T. . 


2yAt\ 

yAt / \ 

f I*!* 1 T’ 1 

AC / 


YAt 

Ir - 

- .) 
1/3-1/ 


2CAC 

\ i/j+1 


yielding 

the first 

familiar 

stability criterion 


At < AC^/2Y 


(8) 


Two modifications of (7) are required for the two boundary condi- 
tions specified in the problem description. 

Insulated boundary . The outer boundary temperatures are de- 
termined from (7) by setting Tj^^j+2 == 'I'i,j-lf thus 



2yAt\ 2yAt 

7?') " 7T 


(9) 


Convection boundary . Using the fictitious node concept posi- 
ted by Ozi^ik [13] , the nodes Ti,^j -1 for j = 1 are defined as 
and their values are arrived at via an energy balance at the con- ^ 
vective boundary 


k 


1 


( 



= h 



(10) 


where Tf^ is the bulk fluid temperature at location ”i". Solving 
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for the fictitious nodal temperatures in (10) yields. 


T, 


f ici 


= 2A^ 


B1(T£, - 




+ T 


i,2 


( 11 ) 


where the Biot number Bi = h R/k]^ Inserting (11) into (7) yields 
the relation 



(12) 


From (12) 


arises a more restrictive time stability criterion. 


AC / 1 BiAC\“^ 

At £ — I — + Bi + ) 

2y \aC / 


to be used in place of (8) , and a second, milder, geometric criter- 
ion from the coefficient of the third term in that 


2 1 

— - — >0 or AC £ 2C-| 
AC Cl 


which implies that the mesh size should be sufficiently smaller 
than the cylinder inner radius. Thus, equations (7), (9) and (12) 

are used to determine all "current" nodal temperatures except T^^q 


3.2. Interface Region 

At this point, Taylor series expansions of temperature are 
imposed from the melt-line to the two known temperatures, Ti^q ^2 
and in the form 


9T 


T, (C) 
k 


= T. 


+ (C - 


e . 

1 


Ef) 


9C 


(C - £.)== a^T 
1 k 


9C 


c . 

1 


k = Ir 2 


(14) 
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generating the following two equations: 


1 ,q-l 


(AC + 

Tin - (AC + 6C.)A. + ^ 

^ ^ 2 




( 15 ) 


'^i,q+2 = + (2AC - 6C.)B. + 


(2AC - 6C^)^ 3^T2 




e . 
1 


(16) 


containing the unknowns Aj^, and the two second derivatives. In 
addition, (17) is valid at the melt-line for all times. 


dTk 

3Tk 

3Tk 

de . 




= 0 = 

+ 

J. 

k = 1, 2 

(17) 

dx 

e . 3x 

1 

e . 3C 

1 

e . dx 
1 




where the first quantity on the right hand side is merely the heat 
conduction equation (4). Thus, (17) gives rise to two more equa- 
tions 


3^T, A. de. 

+ — + A. — - = 0 

3C^ s. ^ dx 


(18a) 


3 




B. 


1 



dx 


0 


(18b) 


with the introduction of another unknown dcj^/dx , which is the Ste- 
fan Condition (5) rewritten below as. 


de 


i 


dx 


1 

A 


(A^ - KB^) 


(19) 


(16) for the second derivative terms 
(18) yields two coupled equations with 


Hence, by solving (15) and 
and substituting (19) into 
the two unknown gradients , 
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0 


( 20 ) 


A. 

1 


B. 

1 


2 


+ 


+ 



. B.K + (t. - Tm^ = 


A.B. 2yX 
-■ — - + 


K K 


^i,q+2) 


0 


( 21 ) 


where = A? + 6^1/ <^2 = 2AC - These two equations are then 

applied at each "ith" location, where A^ and B^ are found via an 
iterative process, such as a Newton-Raphson routine. 


3.3. Determination of and Tj[^q+]^ 

With the known temperature gradients, (15) and (16) can be 
solved for the two second derivatives, in which case everything 
is known about the Taylor-series expansions made near the inter- 
face. Of further significance is the fact that such expansions 
imply a quadratic temperature distribution throughout the inter- 
face regions which is evident in the similarity between (14) and 


T(5) = a + b^ + (22) 


where a = Tm/ b = Aj^ or B^ and c = | . or 1 ^ . 

depending upon the region. Thus, temperatures at^nodes "q" and ^ 
"q + 1" can be calculated from (22) by using C = ~ and ^ = A^ 

- respectively. 

Lastly, the new melt- front velocity at each "i" location is 
computed from (19) . In general, the future melt-line position 
within each pipe section is determined by multiplying the melt-line 
velocity times the time-discretization parameter At (this repre- 
sents an inherent limitation in the finite-difference method where 
the values calculated are to remain constant over the time interval 
At). However, due to the cylindrical nature of the problem, i.e. 
the increase in surface area faced by a radially expanding inter- 
face, the Stefan Condition which heretofore was appropriately 
utilized in the mathematical derivations at the melt-front, must 
be interpreted differently when dealing with a time-increment ana- 
lysis. The energy balance at the interface, assuming unit depth, 
is 


8T, 3T„ pLV 

k, (2Trr. ) — - - k.(2iTr. ) — ^ = (23) 

^ 9r 3r At 


where V is the volume of material undergoing a change of phase dur- 
ing the time increment, or Tr(r|.^^^ - r|.) , and r^ is the radial 
melt-line location at time "t". In non-dimensional form, (23) re- 
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and are useful in providing starting values for the iterative method 
involving successive and Bi calculations. 

Following this starting procedure and while the melt-line is 
within the first two nodes/ the method outlined in sections 3.1 - 
3.3 is used with one modification. Since a value for Ti,q-i is 
non-existent, the Taylor-series expansion is made to node q, over ^ 
the distance -65i, utilizing the temperature Ti,q which at this 
point has not been calculated. An expression for this temperature 
is arrived at by differentiating (14) with respect to the space 
variable and substituting this into the equation resulting from 
an energy balance at the wall (10) , hence 










= Bi 




(28) 


After solving (28) for Ti^l, an identical process to that in 3.2 is 
carried out, resulting in a modified equation to be used in lieu of 
(20) as the melt-line travels between the first two nodes. 




X 

+ 

3Bi 



A.B.K 





(29) 


where B = <SC|/2 + 6^^/Bi 

For the case when the medi\am is initially at the melting tem- 
perature, (Bi = 0) , the starting value for Ai is found from the 
energy balance at the interface, namely 



h (Tjjj 


Tfi) 


or 


hR 

Ai I = -— (Tm - Tf .) = Bi(Tm - Tf .) 
|T=0 


(30) 


and the analysis 3.1 - 3.3 is resumed. 


3.5 Convecting Fluid Analysis 

Figure 4 depicts a horizontal view of the heat exchanger and 
surrounding PCM, showing the smooth-line approximation to the melt- 
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Figure 4. Cross-sectional view of the heat-exchanger with a smooth- 
line approximation to the "M" melt-line locations, and 
an elemental fluid flow section, at an intermediate time 

It I' 

i • 


line locations of the "M" axial sections, as well as a representa- 
tive pipe flow section. An iterative method is employed to deter- 
mine the average bulk temperature Tf ^ , whereby successive energy 
balances are performed on the volume of fluid passing through the 
particular node "i". During the first pass in the method, Tf- is 
assumed to be equal to the incoming fluid temperature Tf^ Thus, 

an energy balance on the fluid element (31) yields the maximum fluid 
outlet temperature Tf^ 

0 =(hA - Tf.) = ™Cj(Tf. (31) 


Of course, the actual bulk fluid temperature Tfj^ is higher than 
that of -j^due to the heating of the fluid as it passes through 

the node, and is now approximated by averaging the newly computed 
values as 




The energy balance (31) is again applied to yield a smaller value 
of and the above process is iterated until successive cal- 

culations of Tf^ are within a desired tolerance. At this point 
Tfi^l . = Tf . and the method is carried out repeatedly until 
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the fluid exits the pipe. The fluid temperatures are adjusted after 
each time step used in the finite-difference routine. 


3.6 Melt-Line Crossing a Node 

When a newly computed melt-line at an "ith" location indicates 
that the melt-line has crossed over node "q+1", an error arises in 
the computation of Ti^q_i (formerly whereby the tem- 

perature Ti,q Ix-At “ '^m+ is erroneously used in the finite-differ- 
ence routine to compute Ti,q-i (Actually, during the time- incre- 
ment, its temperature varies from to Tjn , or, a value less than 

the melting temperature) . To avoid this inconsistency, the time 
step is adjusted (reduced) so that the melt-line falls exactly on 
the node q+1. At this time, the Taylor-series expansions are prop- 
erly carried out, the new melt-line time location is determined by 
resuming the full time step, and the normal-procedure is followed. 


3 . 7 Terminating Procedure 

In the case of an insulated boundary at the following 

routine is used as the melt-line just crosses the node Cn- 1* Simi- 
lar to the starting solution (3.4), the nodal temperature Ti^q +2 is 
non-existent, and the Taylor-series (14) is made only to node ^n« 
Using the insulated boundary condition 


9T, 

— £ = 0 

d 


in (14) where d = A? - solving for the second derivation term, 

and substituting this into ( 18 b) yields 



(32) 


which replaces (21) when the melt-line is between the nodes ?|^_i 
and The current value of T^^N is then computed from (22) which 

reduces to 


T. „ = T„ + 
i,N m 


dB. 


(33) 


3.8 A Third Stability Criterion 

The third and most restrictive stability criterion results from 

an analysis in Appendix I, (adapted from [16]), whereby the start- 
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ing solution used in Sec. 3.4 to generate nodal temperatures is 
meshed with the finite-difference routine yielding the same nodal 
temperatures. The resulting criterion is that the thermal parame- 
ters and node size should be chosen such that B > 1.0 or 


h 

— Ar > 1 


(34) 


Violation of (34) results in starting solutions that yield fluctu- 
ating temperatures and melt-line velocities. 

This condition of course will not arise from the use of a con- 
stant wall temperature at the inner radius. The problem solution 
with this boundary condition can be attained by the method outlined 
in Appendix II, or by the use of convection with a very high value 
for the heat-transfer coefficient, in which case (34) would again 
have to be satisfied and result in an extremely small time-step. 


4 . RESULTS AND CONCLUSIONS 

The following representative program outputs were generated 
from the thermal property data of two types of phase-^change mater- 
ials currently under investigation, namely, salts and metal alloys. 
The purpose of this is twofold: first, to provide usable data con- 

cerning each PCM, and second, to emphasize the ability of the pres- 
ent program in contrasting several PCM's to determine their applica- 
bility. A Burroughs Corporation B-7700 series computer was util- 
ized to illustrate the behavior of three chosen PCM's, the thermal 
properties Of which appear in Table 1. The heat-exchanger system 
is similar in function to that used in [15] where a workable config- 
uration was developed for use in a Solar Brayton Power Plant con- 
cept, and its dimensions appear in Table 2, where an effort was 
made to contrast the materials by choosing equal governing parame- 
ters (e.g. constant h and constant temperature difference between 
the melting and inlet convecting temperatures) . The assumptions used 
in producing the results are listed in Table 3, where only 1 through 
4 are required by the foregoing analysis. 

Figures 5-7 depict the melt-line locations within the heat- 
exchanger at different times, for the salt and metal-alloy PCM's. 
Actually, each interface location represents a contour surrounding 
the internally-convected pipe, so that in effect the three-dimen- 
sional shape of the melt-line resembles a cone in Figure 5, or a 
trumpet end in Figures 6 & 7. Due to the higher conductivity of 
the alloys, more heat is convected within the upstream pipe sections 
than that by the salt, hence, the contours are more curved. The 
slight curving of the fusion fronts with time in Figure 5 is attri- 
butable to the decreasing fluid temperatures within the downstream 
pipe section; as less heat is conducted into the upstream sections, 
its temperature rise is smaller. The lower thermal conductivity of 
the salt yields an expected increase in the time required for solid- 
ification. The similarity of the metal alloys is reflected in the 
similar interface shapes, with the Mg-Si alloy requiring more time 
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Latent 

Specific 


Heat 

Heat 

PCM 

(kj/kg.) 

(kjykg/»C) 

Salt 

598.0 

• 

CO 

o 

7CaF2 - 
54KF - 


i 

39NaF 






Density 

(kg/m.^) 

Melting 

Temperature 

(°C) 

2,339.0 

640.0 

2,250.0 

579.0 

2,000.0 

947.0 


Table 1. Thermal property data for three representative phase-change 
materials . 


H/X inner radius (cm.) 

H/X outer radius (cm.) 

H/X length (cm.) 

Working fluid 

Working fluid specific heat (kj./kg.-®C) 
Mass flow rate (gr./sec.) 

Heat transfer coefficient (w/m.^-®C) 
Initial superheat (®C) 

Working fluid temperature (°C) 


- 2.5 

- 8.0 
- 20.0 

- Helium 

- 5.2 

- 15.0 

- 5454 

- Tm + 


5 

lO'-C 


Tm - 140°C 


Table 2. Heat exchanger data. 


1. Neglect natural convection 

2. Constant melting temperature 

3. Density is a constant in both phases 

4. Neglect axial conduction within PCM and working fluid 

5. Constant material properties 

6 . PCM initially superheated 

7. No crystallization at interface 

8. Heat convection at inner radius 

9 . Insulated outer radius 

10. Linear bulk- temperature distribution within each pipe section 

11. Neglect pipe thermal capacity 

12. Constant heat-transfer coefficient and inlet fluid temperature 


Table 3. Assumptions used to generate results in Figures 5 - 12 . 


to solidify due to its higher latent heat of fusion. Due to the 
small time step (high h value required from (34)) the melt-line 
propagation only up to about 5 cm. was plotted to facilitate compu- 
ter generation of the results; however with greater access to com- 
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5.50 


M 

H 

Ui 


PCM is SALT 


4.75 h 



PIPE LENGTH (cm) 


riME (min) 

a - 0.41 
b - 1.41 
c - 3.09 
d - 5.53 
e “ 8.84 
f- 13.10 


Figure 5. Melt-line location at six different times within the 
heat exchanger for the salt PCM. (Biot No. = 2537.0) 



0.00 4.00 8.00 12.00 16.00 20.00 

PIPE LENGTH (cm) 


Figure 6. Melt-line location at seven different times within 
the heat- exchanger for the metal alloy Al-Si. 

(Biot No. = 145,5) 



INTERFACE . . 

POSITION 



PIPE LENGTH (cm) 

Figure 7. Melt-line location at seven different times within 
the heat- exchanger for the metal alloy Mg-Si. 

(Biot No. = 145.5) 


INTERFACE 

POSITION 



TIME (min) 

a 

- 3.10 

b 

- 8.79 

c 

“ 17.47 

d 

- 29.47 

e 

- 45.50 

f* 

- 67.00 

9 

-95.58 

h 

- 133.12 

• 

j 

- 177.62 


PIPE LENGTH (cm) 

Figure 8. Melt-line location within the heat-exchanger for times up until 
the salt PCM is almost totally solidified. The heat- transfer 
coefficient and mass flowrate v;ere reduced to 731 w/m2-®c and 
1.4 gr/sec, respectively, in order to exhibit the melt line 
propagation along the insulated outer boundary. (Biot No. = 340 
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RADIAL DISTANCE (cm) 


Figure 10. Radial temperature distribution within the first pipe 
section for different times, using the salt PCi^. 

Times are the same as those in Figure 5. 





PCM is Al-SI 


3.60 


RADIAL DISTANCE (cm) 

Radial temperature distribution within the first pipe 
section for different times, using the metal alloy Al-Si, 
Times are the same as those in Figure 6. 


TEMPERATURE rO 



RADIAL DISTANCE (cm) 

Figure 12. Radial temperature distributions at six locations 

throughout the heat-exchanger pipe at an intermediate 
time of 37 seconds, for the metal alloy Al-Si- 


puter process time, the entire solidification region is similarly 
produced. As an example of this. Figure 8 illustrates almost the 
entire region of the salt PCM being solidified, via the reduction 
of both the heat- transfer coefficient and mass-flowrate. The inter- 
face movement along the insulated boundary is shown, as well as the 
increase in solidification time. Also of interest, is the effect 
of mass flowrate on the interface contours. Figures 5 and 8 indi- 
cate that the mass flowrate can be reduced to the point where even 
the salt, with its low thermal conductivity, can transfer much of 
its energy within the upstream nodes, in a manner similar to that 
of the metal alloys. 

The efficacy of the alloys over the salt in charging a heat 
storage system is shown in Figure 9 where the area under the curves 
represents the total energy transferred after a certain time, and 
is admeasure of the effectiveness of the heat exchanger. The values 
for Q are determined by the mass flow rate, m, specific heat of the 
fluid, Cf, and the temperature difference between the inlet and 
outlet fluid temperatures. Of course Q max occurs when the tempera- 
ture difference is the greatest (i.e. initial temperature of the 
medium minus the inlet fluid temperature) . 

Figures 10 and 11 exhibit the radial temperatures within the 
first pipe section at various times, for the salt and Al-Si data 
respectively. The dissipation of heat with time throughout the 
superheated region is evident, as well as the temperature non-lin- 
earity at the interface. The convecting surface temperature is also 
depicted as decreasing with time. The fact that the inside surface 
temperature of the pipe is lower for the salt than the alloy, and 
not vice versa as one would be lead to believe, is due to the heat 
convection value chosen for the analysis. Similar to the resistance- 
network concept, the convection resistance at the wall is greater 
than that encountered within the alloy medium, but of the same order 
of magnitude as that within the salt. Hence, the rate-determining 
step for heat transfer out of the metal alloy, is the convection 
resistance (unlike the salt) , and if this were made smaller (by in- 
creasing h) then the surface temperature would decrease. In a 
sense, the heat-transfer efficiency of the alloys is even greater 
than that shown. Figure 12 provides an axial composite of radial 
temperature distributions throughout the partially solidified metal 
alloy Al-Si at a time of 37 seconds. The inlet and outlet pipe 
sections are 1 and 11, and every other section is included. The up- 
stream nodes are undergoing solidification, while the downstream 
sections are approaching the melting temperature. 
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NOMENCLATURE- (non-dimensionalized parameters appear in parenthesis) 

temperature gradient at interface in solid region 
temperature gradient at interface in liquid region 
Bi Biot number, h R/kx 

Cx specific heat of solid phase 

specific heat of working fluid 
h heat- transfer coefficient 

kx thermal conductivity, solid phase 

k 2 thermal conductivity, liquid region 

K ratio of kl/k2 

L latent heat of PCM 

a axial length of pipe section 

m mass flowrate 

M number of pipe sections 

N number of radial nodes 

q node in solid region closest to interface 

q+1 node in liquid region closest to interface 

R- reference length 

Rl(^l) pipe inner radius 

1^2 (^N> outer radius 

r(^) radial coordinate 

Rj(Cj) radial location of node 

distance between node q and melt-line 
Ar(A^) nodal spacing, R 2 “Rj^/N-l 
Si(£i) radial distance to interface 
t(T) time 

At (A t) time increment 

T_. fictitious node 

fxcx 

T, fluid temperature 

T fusion temperature 

wall temperature 

T. . nodal temperature; l<i<M;l<j<N 
thermal diffusivity, solid region 
thermal diffusivity, liquid region 
p density of PCM 

A non-dimensionalized parameter, L/C^ 

Y equals within liquid region; equals 1 within solid region 
y^ ratio of thermal dif f usivities 
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APPENDIX £. determination OF STABILITY CRITERION 

An energy-balance at the convection boundary during the onset 
of solidification yields. 


3T 

= Bi(T, - T.) ( 1 ) 


where the subscript "i" has been dropped for simplicity, and (25) 
provides an expression for 3T/85 |r • Solving for the temper- 
ature at the boundary yields i-O 


4T, - T, + 2BiA? T^ 


3 + 2BiAC 


Now after an infinitesimal time step At, (2) should also hold for 
the new values of Tj^, T 2 and T 3 , denoted by primes as, 

, 4t' - t' + 2BiA? T, 

T = — f e £ (3) 

3 + 2BiA^ 


t I 

but T 2 and T 3 are found from the finite-difference equation 


I 


T 2 = T 2 

a + b (T^ + 

Tf) 

+ c(T3 - 

Ti) 

1 

T = T 
^3 3 

a + b(T^ + 

T 2 ) 

+ c(T^ - 

T 2 ) 

where a = 1 - 

■ 2yAt/A5^ 

b = 

YAt/A5^ , 

c = Aty/Z^AC 


Svibstituting (4a, b) into (3) yields 


(7) , as 
(4a) 

(4b) 


T^ = l/d[Tj^(4b - 4c) + T2(4a - b + c) - T^ (a - 4b - 4c) 

+ T^(-b - c) + 2ACBiT^] (5) 

where d = 2A5Bi + 3, and 
(5) can be written in simpler form as 

Tji^ = l/d[oTj|^ + 3 T 2 - bT 3 - (j)T^ + co] (6) 
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where 



9yAT Aty 

e = 4 + 

6yAt 4Aty 

9 = 1 ; 

2CAC 

YAt Aty 

= — T + 

AC 2CAC 

(U = 2ACBiT^ 


(7a) 

(7b) 

(7c) 

(7d) 

(7e) 


The error obtained in (6) is a result of individual errors in all 
the computed temperatures, hence 


AT^ = [ctAT^ + 3 AT 2 - nAT^ - 4>AT^]/d 


(8) 


and the maximum error obtainable from (8) occurs when all the tem- 
perature errors have the same sign, hence. 


at' = = [|a||AT I + iBllAT-l + InllAT^I + |(|)1 IAT 4 I ]/ d 

(9) 


Since the temperature Ti changes the most due to the heat convec- 
tion at the boundary, its error should be the largest of all asso- 
ciated errors, hence. 


at' = At 1 «i/d [|a| + |8l + |n| + 
^ max ^ ^ 


( 10 ) 


I 

If 6 , in (10) is less than ^ the right hand side of the equation, then 
the successive values of T^^ at later times is assured to contain 
an error converging to zero, and in a sense a bound is determined 
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for the maximum tolerable error used in the finite-difference equa- 
tions. From (7d) (f) is always positive, and (7a) contains the second 
stability criterion where A? « so a is always positive. Thus 

(10) alters to 


^1 — '^ i /^ 1^1 


I 

and stability is achieved when 6^^ < which implies 
o + <{)+ |3| + |n| <d 


( 11 ) 


(12) 


Rewriting (7) with v = Ary/A^^ yields 


AC 

a = 4v[l - — ] 

2C 


vAC 

3 = 4 - 9v + 

2C 


2ACv 

n = 1 - 6v - 

K 


vAC 

(j) = V + 

2C 


(13a) 


(13b) 


(13c) 


(13d) 


Case 1 . II is positive, which implies v(6 + 2AC/C) < 1 or v < 
1/6 in which case then 3 is surely positive and (12) is then sim- 
plified to 


1 - 5v < ACBi or ACBi > 1 (14) 


which would satisfy the above criteria. 

Case 2 . 3 is positive, which implies v < 4/9, hence the abso- 
lute value of n would be 6v + 2ACv/C - !• After plugging into (12) 
and simplifying, this leads to v < ACBi which is satisfied also 
by BiAC > 1. 
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APPENDIX II. CONSTANT WALL TEMPERATURE ANALYSIS 


The constant-wall temperature boundary condition represents 
a more complicated and paradoxical condition for use in phase- 
change problem initiation where the medium is initially at a tem- 
perature different from the fusion temperature (i.e. superheated) . 
Essentially, this condition states that at time equal to zero, the 
wall temperature is at a temperature below the fusion point. This 
ignores that fact that the wall temperature is actually at the 
initial temperature at t = 0, and also that a finite amount of time 
is required to reduce it to the melting temperature, while also 
allowing a non-constant temperature distribution to establish with- 
in the superheated region. Thus, any attempt to model such a pro- 
cess is at best, an approximation to an unrealistic condition. 

One approach is to choose an initial time step significantly 
smaller than the one computed front Equation 8, Section 3.1, and 
perform an energy-balance on the small volume of material undergo- 
ing phase- change during that time increment. 




A<T„ - T„) 


6r 


\ At = pA«r(L + - v) 


(i) 


and the non-dimensional distance 6? is solved from (1) as 


6 ? = 


(5r 


R 


^ / ^1 ^'^m ~ ^w^ \ 


1/2 


where the cylindrical volume effects are considered negligible for 
such a thin ring of material. Following the establishment of an 
initial 6^, the remaining procedure is similar to that of sections 
3.1 to 3.3, with the simplification that Tj^ in Equation (15), 

Section 3 .2 for i = 2 (or Tj^ „ when the interface is between the 
first two nodes) is known as rhe constant T^^^. Lastly, had the ini- 
tial temperature distribution been something other than a constant, 
such as a quadratic or parabolic, the starting solution outlined 
in 3' . 4 would apply just as well, and obviate the inherent diffi- 
culty in determining the solid temperature gradients, A^, first. 
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INTRODUCTION 


Materials which undergo a change of phase are being intensely 
researched for applications in the fields of Thermal Energy 
Storage and Transfer. The high energy-density and desirable 
thermal characteristics of these phase-change materials (PCM's) 
make them attractive candidates for use in such areas as 
solar energy collection and peak- load leveling [20] . Along with 
this increased interest in PCM's is the need for accurate 
modeling of the thermal behavior of such materials. Due to the 
inherent temperature nonlinearity at the melt-line, the analyt- 
ical solutions available for solving the phase-change process 
are relatively few in number, and require restrictive boundary 
or initial conditions [1-8] . The numerical approaches 
alternately employed for the problem solution, although some- 
what more applicable, also require restrictive conditions which 
many times reduces the analysis to that of a pseudophase-change 
process far removed from the actual system itself. Extehsion 
to two-dimensional cyclindrical coordinates further reduces the 
number of available solutions in the literature [9-15] . Talwar 
and Dilpare [9] appears to be the first attempt at modeling 
such a phase-change problem. 

The method presented in this paper evolved from that initially 
posited by Forster [10] , refined by Murray [11] , and subse- 
quently utilized by Lazaridis [12] , resulting in a complicated 
routine difficult to implement in cylindrical coordinates. The 
feasibility of this solution was shown by Farouk [19] , but the 
required clarification and inherent simplicity was demonstrated 
by Harding [16] and Labdon [17] . 

The proposed solution method represents an efficient numerical 
routine for modeling two-dimensional heat conduction (axial 
and radial) within an annulus of phase-change material. 
Finite-difference equations monitor the heat conduction in the 
regions away from the interface, while truncated Taylor-series 
expansions of temperature are made near the interface, so as 
to couple the thermal behavior between the solid and liquid 
regions. In this manner, the exact melt-line locations and 
temperature histories are continually determined as the 
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interface propagates throughout the medium. Inputs include time 
and spatial-dependent boundary conditions, temperature-dependent 
properties and initial temperature distributions, while outputs 
include interface positions and temperature distributions as 
functions of time. 

PROBLEM STATEMENT 

The cylindrical geometry 
under consideration is 
depicted in Figures 1 
and 2. Figure 1 is an 
axial view of a hollow 
cylinder of length ' , 
containing a PCM parti- 
tioned into 'N' equal- 
thickness, concentric, 
annular regions or nodes, 
showing a particular 
melt-line location. For 
the following analysis, 
heat convection is 
imposed on the inner 
cylinder R^, while two- 
dimensional solidifica- 
tion occurs radially 
outward toward the 
insulated boundary R 2 . 


The concatenation of 'M' cylinder sections results in a heat 
exchanger of desired length 'd', as illustrated in Figure 2. An 
energy balance is performed upon the bulk-temperature of the 


PCM 



Figure 2. Horizontal view of the hollow cylinder of length 'd', 
where d = Mi,, containing a PCM between an internally convected 
pipe and an insulated outer one. 



Figure 1. Axial view of the cylinder 
with representative concentric rings 
(nodes) of PCM, depicting a particular 
melt-line location. 
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convecting fluid as it travels through each cylinder section. In 
this manner, the melt-line contour is continually traced through- 
out the phase-change medium, as illustrated in Figure 3, which 
also shows the two-dimensional mesh network. 

PROBLEM SOLUTION 

The basic assumptions for the analysis are that the fusion temper- 
ature is constant, no crystallization occurs at the interface, 
the density is constant in both phases and that natural convection 
effects are neglected. The governing equations and conditions used 
throughout the medium are : 



Figure 3. Cross-sectional view of the cylinder with a smooth- 
line approximation to the 'M+1', radial melt-line locations, and 
an elemental fluid flow section, at an intermediate time 't'. 
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3T 

k 

II 

1 

II 

= 0 

(4) 

3r 

R 

N 

0 

II 

N 

z=d 


T = 

T r=S(t) 

m 


(5) 

T = 

(z,r) 

t = 0 


(6) 


where k = s or Jl depending upon the solid or liquid region. The 
following non-dimensionalizing parameters are introduced into 
the above equation. 


C s r/R ; T s Og t/R^ ; e(x) 5 S(t)/R ; A? s Ar/R ; K 5 k^^/k^ 
^ E R^ ; Cjj = R2 C = z/R 


(where R is a chosen reference length) generating the following 
equation in which all variables but temperature are non- 
dimensionalized : 


3^T, 1 3T, 1 3T, 

k ^ k ^ k k 

9C^ C 3^ Y 3t 


1 i < e, k = s 
^ ^ k = Z 



5 = e (t) 


(7a, b) 


(8) 


3T 


H 


= Bi(VTf) 


5 = Cl / T > 0 (9) 


_ 

= 'Zk 

1 = 0 

0 

A 

(10) 


3C 

II 



T„ 5 = e(T) 

; T, = 

T^(C,C) 

X - 0 

(11,12) 
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Method 1 - One-Dimensional Taylor Series Expansion 


Figure 4 describes the temperature distribution within the inter- 
face region for a particular cylinder section. The axial temper- 
ature distributions in the + c direction (not shown) account 
for the two-dimensional movement of the interface. At this point 
all that is known are the 'current' melt^line locations at time t 
and the 'old' nodal temperatures at time t-At. The object of the 
method is to determine the current nodal temperatures and new 
melt-line velocity, which requires the successful execution of the 
following computational logic ; 

Step 1 - Determine Melt-Line Slopes . The more arduous task 
of monitoring fusion- front propagation along two orthogonal 
directions is superseded by the method of determining the radial 
projection of the normal velocity gradient, shown graphically in 
Figure 5. The mathematical derivation of this procedure appears 

in [13], the result of which is a factor |l + {^) = to be 

multiplied by the normal gradient, that allows contznuous moni- 
toring of the melt- line gradient in the radial direction. The 

3 s 

term -r— represents the slope of the melt-line in the C direction, 
and is determined at each 'i' location as / (2 A^ ) . 


T 



Figure 4. Radial temperature distribution near the interface, 
within a particular pipe section. 



e 



Figure 5. Particular melt-line location depicting radial projec- 
tion of normal directional gradient, and distances used for 
localized melt-line slope determination. 


Step 2 - Nodal Temperatures . The current nodal temperatures 
, j where l_<ij<M; lf.;jjf.q“l and q + 2 £ j £ N are computed 
via*^the finite-difference equation 


1 . = T. ■ 1 

i/D i/D \ 

2yAt 2yAt^ 

( AC^ AC=" ! 


+ YAt I 
2CAC ' 

^^i,j+l " '^i,j-lj 



with an associated time stability criterion, 


At 


< /_i_ + _i_ - Bi Bi y ^ 

2y 1 2K AC j 


(14) 


determined from the convection boundary [18] 
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step 3 - Determine Axial Temperatures and Gradients ♦ The 
four axial, inter-node temperatures shown in Figure 6 (dashed 
circles) are required (at each location) for determining 
temperature gradients and second-derivatives at the melt-line. 
Each temperature is interpolated from a three-point distribu- 
tion determined from the radial temperature profiles originating 
from the interface (solid circles) . Using these inter-node 
temperatures and the melting temperature, two quadratic profiles 
in the form 


T (c) = a + be + cC^ 


(15) 


are utilized, both of which originate at the interface and extend 

in the ± e directions (horizontal dashed ellipsis of Figure 6) . 

The two second derivatives of temperature in the axial direction 

and then determined at the melt-line from (15) and are referred 

to as Z " , Z " . 
s' I 


Step 4 - Determine Radial Temperature Gradients . Radial 

Taylor-series expansions of temperature are then applied from 

each melt-line location in the ± ? directions (vertical 

dashed ellipsis of Figure 6) to the known value of nodal 

temperature T. ^ , and T. in the form: 

i,q-i i,q+2 


3T, 


Tj^(c) = Tj^ 


"i 


+ (? - e^) 




(5 - £i) 


^i 


^ 9^T, 


3^ 


(16a, b) 


"i 




Figure 6. Local tem- 
peratures and distri- 
butions required for 
temperature derivative 
calculations at a 
particular melt- line 
location. 
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In addition, (17) applies at the interface 


dx 




3 


(17a, b) 


The seven equations (7a, b) , (8), (16a, b) , (17a, b) are then 

combined in a manner similar to that of [18J to solve for the 
seven unknowns, namely , the two second derivatives of temperature, 
the two radial temperature gradients, the two time derivatives 
of temperature, and the melt-line velocity, all of which are 
applied at the current interface location. 

9T 

equations are coupled in the unknowns — 

A. and B. , 

X X 


The resulting two 
3T, . 

defined as 


and — 


■£ 





step 5 - Calculate Interface Region Nodal Temperatures . The 
similarity between (15) and (16) indicates that the Taylor-series 
expansion of three terms is identical to instituting a quadratic 
temperature profile, and thus the foregoing analysis of the inter- 
face regions assures uniformity throughout the method. At this 
point, (16) can be used to determine the remaining nodal termper- 

atures T . and T . . , . 

x,q x,q+l 

Step 6 - Determine New Melt-Front Velocity . Via an analysis 
in [18] the Stefan condition in cylindrical coordinates (8), which 
heretofore was appropriately utilized in the mathematical deriva- 
tions, is interpreted differently for use in a numerical routine. 
The increase in surface area faced by a radially expanding inter- 
face alters (8) by the factor below 

. * + t \ 

T + AT T_ 1 

/ 



A. 

X 


KB. 

X 


( 20 ) 



which reduces to unity as the time-discretization parameter 
At is made smaller. 


Method 2 - Two-Dimensional Taylor-Series Expansion 


As shown in Figure 7, 2 two-dimensional Taylor-series expansions 
of temperature are applied at the interface in the ± ^ and ± ? 
directions in the form: 


9T 

T(c,C) = T(i,e.) + — 


(C-i) 







+ 


8^T 

35 ^ 


^i 




2 


3^T 

+ 


(5-e. ) (c-i) 
e. 

X 


= + a(c-i) + b(C-£i) + cCc-i)"" + d(?-e^f + e(c-i) (C-e^) (21) 


where the constants represent the appropriate gradient or second 
derivative term in (21) / and the terms in parentheses represent 


€ 



Figure 7. Temperature field surfaces generated from two-dimensional 
Taylor-series expansions about the interface. 


F9 



the axial and radial distance to a particular node of known 
current temperature. Thus, the selection of five known temper- 
atures in the neighborhood of the interface, in both the liquid 
and solid regions, is sufficient to determine the sought after 
temperature gradients and in (21) . Intuitively, the 
closer the chosen nodal points to the interface, the more accuracy 
that would result in generating the two-dimensional surfaces of 
temperature, and an intricate computer program was utilized to 
determine the nodes of closest proximity to the interface. Thus, 
steps 3, 4 and 5 of Method 1 are superseded by the above proce- 
dure, which results in decreased computer process times for 
problem solutions. 

Information as to starting and terminating procedures, the 
fluid energy balance, the melt-line crossing a nodal point, and 
the use of a constant wall temperature boundary condition can be 
found in [18] . In particular, the starting procedure requires 
the following restrictive stability criterion. 



Ar > 1.0 


RESULTS 


The following system parameters were chosen to describe the 
various outputs from the method: 


Cylinder inner radius (cm) 

Cylinder outer radius (cm) 

Cylinder length (cm) 

Working fluid 

Working fluid specific heat (kj/kg-“C) 
Mass flow rate (gr/sec) 

Heat transfer coefficient (w/m^’-’C) 
Initial superheat (*C) 

Working fluid temperature (®C) 


- 2.5 

- 8.0 

- 20.0 

- Helium 


- 5.2 

- 15.0 

- 5454.5 


- T„ + 10 

- T? - 140 

m 


In addition, the following assumptions were used to facilitate 
data generation only, and are not required by the foregoing 
analysis : 


1 . 

2 . 

3. 

4. 

5. 


6 . 

7. 


Constant material properties 
PCM initially superheated 
Heat convection at Inner radius 
Insulated outer radius 

Linear bulk- temperature distribution within each pipe 
section 

Neglect pipe thermal capacity 

Constant heat-transfer coefficient and inlet fluid 
temperature 
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Figures 8 and 9 depict melt-line contours at various times for 
the salt PCM 7CaF2“54KF-39NaF and the metal alloy PCM Al-Si, 
respectively. The former, having a lower thermal conductivity, 
yields the result of relatively straight contours and longer 
times required for solidification. Figure 9 illustrates the 
more curved interface shapes due to the greater heat-transfer 
within the upstream nodes. 

Figure 10 contains radial temperatures within the first cylinder 
section at times corresponding to those in Figure 8 for the salt 
PCM. The decreasing inner surface temperature and initial 
superheat with time are evident. 

Figure 11 represents an axial composite of several radial 
temperature distributions within the Al-Si PCM at a time of 
20 seconds following the onset of solidification. The upstream 
sections are undergoing phase-change while the initial super- 
heat is still being dissipated within the downstream ones. 
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Figure 8 - Melt-line contours at six different t 
the salt phase-change material 7CaF2-5 
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Figure 9 - Melt-line contours at seven different times within the annulus 
for the metal alloy phase-change material Al-Si. 
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Figure 10 — Radial temperature distributions within the first cylinder section 
for the salt PCM, at times corresponding to those of Figure 8. 
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Figure 11 - Axial composite of radial temperature distributions within the 
cylinder at a time of 20 seconds following the onset of 
solidification of the metal alloy PCM. 


NOMENCLATURE (non-dimensionalized parameters appear in parentheses) 


A. 
1 

B. 
1 

Bi 


C 


1 


C 


f 


h 


K 


L 

I 

m 

M 


N 


q 

q+l 

R 


r(S) 

R. (C • ) 

3 ^3 

6r.(6?.) 
Ar (AO 


S. (e . ) 
1 1 


temperature gradient at interface in solid region 
temperature gradient at interface in liquid region 
Biot number, h R/k 

s 

specific heat of solid phase 
specific heat of working fluid 
heat-transfer coefficient 
thermal conductivity, solid phase 
thermal conductivity, liquid region 
ratio of k„/k 

s 

latent heat of PCM 

axial length of pipe section 

mass flowrate 

number of pipe sections 

number of radial nodes 

node in solid region closest to interface 

node in liquid region -closest to interface ‘ 

reference length 

pipe inner radius 

pipe outer radius 

radial coordinate 

radial location of node 

distance between node q and melt-line 

nodal spacing (R 2 -R-j^) / (N-1) 

radial distance to interface 
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t (t) 
At(Ax) 




z (^ ) 


Z 


It 

k 


time 

time increment 
fluid temperature 
fusion temperature 
wall temperature 

nodal temperature; l£i£M ; l£j 
axial coordinate 

second derivatives of temperature in the axial 
direction evaluated at the interface 



distance from node q to melt-line 
distance from melt-line to node q + 1 


thermal diffusivity, solid region 

0.2 thermal diffusivity, liquid region 

/ 3 s 

B^ = l+(-^l which is projection of normal direction onto 
radial coordinate 

p density of PCM 

X non-dimensionalized parameter, L/C^^ 

Y equals y ^ within liquid region; equals 1 within 

solid region 

y ^ ratio of thermal dif fusivities 02/0.2, 
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CONTAINMENT OF HEAT STORAGE ALLOYS 


Nikhil S. Nagaswami 




ABSTRACT 


The present work is directed towards finding 
optimum containment materials for several recently 
discovered heat storage alloys. These alloys include 
binary, ternary and quarternary systems of the elements 
aluminum, copper, magnesium, zinc, phosphorous and 
silicon, all fabricated or synthesized at the University 
of Delaware by Alan Riechman and Diana Farkas . The 
commercial viability of these systems hinges on their 
proving more cost effective as compared to several molten 
salt systems in their respective temperature ranges. 

Since it is anticipated that the containment system could 
prove to be as or more expensive than the storage media, 
this area of research warrants close study. 

The containment materials must satisfy several 
requirements, amongst them: good thermal stability in 

the temperature range, thermal shock resistance, high 
strength at elevated temperatures, and remain inert to 
attack by the alloy, that is, they must exhibit limited 
solubility, resistance to embrittlement, little or no 
decomposition, resistance to infiltration or liquid metal 
penetration (unless the infiltration enhances strength) and 
form a very slowly growing reaction layer that is protective. 


G1 



Experimental studies completed and in progress 
include: 

A study of wetting, reaction layer growth, and alloy 
penetration and diffusion into the containment material, 
upon exposure, without pressure at temperatures 100 to 
150°C above the melting points of the alloys. 

Similar studies are also being done for the 
individual elements comprising these alloys. 

The experimental techniques employed in the 
present work include optical metallography, scanning 
electron microscopy and electron probe microanalysis. 

Six alloys systems are presently being studied. 

These include three alloys which show good heat storage 
properties on a volume basis and three alloys that are good 
on a weight basis. These alloys show promise for application 
in three temperature ranges: 550-625°, 700-800°, and 
950-1050°C. 

The containment materials being studied include 
two stainless steels, S.S.304 and S.S.316L, and a cheap 
silicon carbide compact, "Carbofrax SiC,” from the 
Carborundum Company. 

In all the systems studied thus far, the SiC shows 
great promise as a containment material, while the stainless 
steels are attacked to a greater degree by the molten alloys. 
A conclusive report on all studies conducted will be 
available soon. 

G2 



INTRODUCTION 


An extensive literature survey has been conducted 
and various aspects of the containment problem studied. 

These include; liquid metal and solid-vapor phase 
embrittlement of metals and refractory structural 
materials, including mechanisms of embrittlement, 
specificity of embrittlement, rates of embrittlement, 
and nature of embrittlement; reaction of liquid metal- 
structural material couples and modes of decomposition, 
attack and material diffusion, surface phenomena and morpho- 
logies and effect on liquid metal interactions, 'infiltration 
of metals and refractory materials by liquid metals, silicon 
carbide — its preparation, properties and reactions, 
stainless steels — their properties and the effect of liquid 
metal environments, analytical procedures and theory 
behind SEM, TEM, STEM and diffraction studies and applica- 
tions to the present problem, heat storage alloy applications 
and future trends for research. 

The conclusions of the literature search are 
presented in the introduction to the thesis based on the 
present work. 
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EXPERIMENTAL DETAILS 


Alloy Preparation 

The first method utilized in preparing the alloys 
was melting in a resistance furnace under argon atmospheres. 
The alloying elements were mixed in a graphite crucible 
which was closed with a tight fitting lid, incorporating 
a stirrer to facilitate mixing of the molten alloy. The 
tight fitting lid was used to minimize vaporization losses, 
particularly of the elements Zn, Mg and p. The set up is 
illustrated in Figure 1. 

The temperature of melting depended on the volatile 
elements Zn, Mg. P present and on the highest melting 
constituent present. The graphite crucible containing 
the elements was introduced into the furnace which is 
preheated to the melting temperature desired and flushed 
with argon prior to and during the melting process. The 
alloy is cooled in the argon atmosphere by switching off 
the furnace and maintaining the argon flow. 

This method proved largely unsuccessful owing to 
several reasons: 

1) When the temperature was maintained at the 
melting temperature of the highest melting constituent, 
excessive vaporization of the volatile elements occurs. 
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2) When the melting temperature is held intermediate 


between the melting temperatures of the volatile and high 
melting elements, incomplete melting occurs and the alloy 
is extremely heterogeneous- 

3) Owing to the air trapped in the crucible 
within the alloying element mix, considerable oxidation of 
the alloying elements occurs, resulting in incomplete 
melting and mixing and in a porous sintered mass, even 
at temperatures higher than the melting temperature of 
the highest melting element. 

In seeking a solution to the above problem the 
following set up was utilized (Figure 2) . 

The alloying elements (in the form of chunks' — 
not powders) were weighed and mixed in a graphite crucible. 
Which is placed within a quartz tube as shown. The quartz 
tube is fitted with an air tight seal and evacuated at 
room temperature to remove almost all the air in the 
alloying element mix. The vacuum tight valve is then closed 
and the tube flushed with argon. In this state, the quartz 
tube placed within a set of water cooled Cu induction 
coils is heated. An argon pressure in the tube is maintained 
by regulating inlet-outlet pressures, while the graphite 
crucible and "chunky" elements provide for good coupling 
with the induction coils. Temperatures are measured by 
means of a thermocouple and an optical pyrometer. 

Following the melting, the alloy is cooled under argon. 
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Fig Z 
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This system provides for quick, convenient melting 
of the alloys and minimizes vaporization losses. In 
addition the induction furnace set up provides good mixing 
of the alloying elements, while minimizing oxidation due 
to trapped air. 

Containment Study 

The experimental set up utilized in this study 
is illustrated in Figure 3. 

A tube furnace constructed by winding Kanthal 
resistance wire around a high grade alumina tube and packed 
with refractory firebrick insulation is used. A high 
density mullite tube is mounted within the alumina tube 
to provide a contamination-free environment (with minimum 
outgassing of the heater tube) and to provide more uniform 
temperatures in the effective reaction zone. A gas inlet 
port is cemented at one end of the mullite tube. At the 
other end a nearly gas tight closure has a sliding face 
shield to facilitate insertion and removal of the specimen 
holder and temperature measurement. Tantalum baffles are 
inserted into the mullite tube between the specimen 
holder and the ends to minimize heat loss to the ends 
and to provide a uniform temperature zone over the length 
of the specimen holder. The construction of the specimen 
holder is illustrated in Figure 4. The entire holder is 
constructed of 304 SS and is free of warpage upon heating and 
cooling. The alloy/base material samples are placed in 
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individual receptables in the holder and are removed by 
means of a tight fitting SS304 strip bent to fit the 
receptacle. Each sample is further isolated in its 
receptacle by means of a tight fitting SS304 sheet which 
slides over the holder, thus preventing excessive 
volatilization and loss of elements and contamination of 
samples . 

The reaction temperature is maintained 100-200°C 
above the melting point of the alloy, to allow enhanced 
reaction/diffusion/penetration rates and ensure complete 
melting of the alloys. Volatilization of the elements from 
the alloy is much lower than it is from the elements 
themselves and does not require using the minimum temperatures 
required for melting. 

The specimen holder with the samples isolated by 
the covering strip is introduced into the furnace which is 
flushed with argon and maintained at temperature - 
The furnace is closed off, and an argon flow sufficient 
to maintain a pure atmosphere is provided. The durations 
of the exposures used are 24, 48, 100, 300 hours, done in 
individual batches. Following each heat the specimens 
are cooled to room temperature in the argon atmosphere by 
switching off the furnace heat. The samples are then 
removed carefully and identified prior to analysis. 
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Analytic Procedures 

The samples obtained in the containment study 
(alloy element or alloys on base material) are visually 
inspected for extent of wetting and reaction. 

In those instances that the alloy wets the 
base or forms a sessile drop on cooling, the effect 
is recorded using an optical metallograph. 

The alloy is then ground off the base material 
to within 0.5mm of the surface, and the base is sectioned 
by means of a diamond impregnated cutoff wheel. 

Different grades of abrasive wheel are used for the SS304, 
316 and SiC samples. The cutoff position is selected 
to obtain a representative section. 

The sample is then mounted in epoxy resin (cured 
at room temperature) . Alumina fillers are used to improve 
the hardness and abrasion resistance of the mold to more 
closely match the sample, thus providing better edge 
retention and minimizing relief between sample and mount. 

The samples are then ground on SiC papers down to 
600 mesh and subsequently polished on diamond impregnated 
cloths (low nap for better retention of inclusions and 
precipitates) for SiC and SS304 and 316. Alumina polishing 
abrasives are used also in the case of the stainless steels. 
The final grade of abrasive used is down to 0.25ym. 
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The sample surface is thus polished flat and free 
of surface relief, well suited for SEM observation and 
EDAX analysis. 

Problems associated with smearing the softer 
alloying elements are thus minimized and practically 
eliminated by grinding off excess alloy, sectioning 
through the unattacked rear face of the sample, grinding 
and polishing while keeping the specimen reaction edge 
suitably oriented, and cleaning efficiently by means 
of swabbing with alcohol, using ultrasonic baths between 
preparation stages and at the end of preparation. 

Etching is only carried out on as-received samples 
of the base materials and on base materials subjected 
to the same thermal cycles under argon as the alloy/base 
material samples, to remove artifacts due to thermal 
cycling alone, such as alloying element segregation. 

SEM micrographs are used to study the growth of 
surface reaction layers and penetration-infiltration of the 
base material, while EDAX-electron microprobe analysis, 
performed on beam-size spots and in-line scans at successive 
depths into the base material from the reaction layer, 
are used to study elemental diffusion and distribution. 
Suitable corrections (ZAF, magnification, specimen tilt 
with respect to electron beam and detector window) are 
applied and spectral artifacts eliminated prior to analysis 
of the results obtained. 
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This system of analysis is seen to provide an 


effective, reliable means of studying the containment, 
problem. 

Preliminary qualitative results are included in 
Figures 5a-f. 

SUMMARY 

The alloy systems- used are listed in Table 1 along 
with related physical and thermal properties . 

The base materials studied include SS304, SS316L 
and Sic ("Carbofrax" from Carborundum Co.). 

Results obtained to date indicate considerable 
reactivity of the stainless steels towards several of the 
heat storage alloys, while the SiC remains inert or forms 
a very slowly growing reaction layer with the molten 
alloys. The surface oxide on SiC has been found to be 
principally responsible for the poor wetting characteristics 
exhibited by the alloys. However; since several examples 
are cited in the literature wherein the liquid metal 
finally diffuses in through a resistant surface layer to 
attack the bulk or undermines the surface layer by initial 
infiltration at a crack, additional studies are being 
undertaken wherein the Si02 layer on the SiC samples is 
removed to the extent possible by means of treatments with HF. 
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Figure 5. Reaction of stainless steels and silicon carbide 
with molten heat storage alloys. Elemental 
diffusion curves are based on the longest 
exposure indicated. Concentrations are relative. 
The apparent effective thickness of the reaction 
layer and its dependence on time is shown at the 
top of each curve. 
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Additional areas of research worthy of study in 
finding a solution to the present problem include: 

a) Special si: 'face or heat treatment of the 
stainless steels. 

b) Studies on SiC surface layer structures and 
the dependence of liquid metal interaction on this aspect. 

c) The development, of SiC coatings and their use 
with stainless steels or other suitable structural material 
to overcome certain commercial drawbacks in the shaping, 

of complex SiC shapes. 

d) Containment studies carried out over long 
periods and under pressures such as that to be found at 
the bottom of a tall liquid alloy bath. 

e) Effect on mechanical properties of the base 
materials due to reaction/infiltration by liquid metals. 
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